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Section  1  Introduction 


The  word  anemometer  stems  from  the  Greek  "anemos"  meaning  "wind" 
and  "metron"  meaning  "measure".  A  hot-wire  anemometer  senses  any 
changes  in  the  variables  which  affect  the  rate  of  heat-transfer  between 
the  wire  and  the  fluid.  Variations  in  heat  transfer  coefficient  can  change 
both  wire  temperature  and  resistance.  If  the  wire  is  made  part  of  a 
suitable  electrical  circuit,  these  changes  can  be  used  to  generate  a  signal 
which  is  related  to  the  instantaneous  heat  transfer.  Thus,  as  Morkovin 
(ref.  1)  points  out,  for  correct  interpretation  of  the  electrical  signal  we 
need  to  know:  1)  the  identity  of  possible  fluid  flow  variations  (eg. 
turbulence  or  sound),  2)  the  laws  of  heat  transfer  between  the  wire  and 
fluid,  3)  the  variation  of  wire  resistance  with  temperature  and  the  effects 
of  conduction  to  the  supports,  and  4)  the  response  of  the  associated 
electrical  system  which  produces  the  measured  current  or  voltage 
variations. 

Unfortunately,  our  knowledge  in  each  of  these  categories  is  far  from 
complete  and  could  well  be  responsible  for  the  current  lack  of  reliable 
data.  A  review  of  hot-wire  data  taken  in  zero  pressure  gradient, 
adiabatic  or  isothermal  wall  boundary  layers  illustrates  the  problem. 
Figure  1.1  shows  data  from  several  sources  for  the  fluctuating  axial 
velocity  component.  The  scatter  is  so  large  tha*'  it  is  impossible  to 
construe  that  any  form  of  similarity  with  Reynolds  or  Mach  number 
exists.  The  picture  is  even  more  confusing  when  the  distributions  of  the 
other  two  normal  stresses  are  plotted.  Fig.  1.2  shows  the  measured  shear 


stress  distributions  and  again  no  pattern  of  similarity  can  be  observed. 
Indeed,  only  Klebanoff's  incompressible  measurements  (shown  for 
comparison)  approach  the  anticipated  limiting  value  of  unity  in  the  wall 
region.  These  results  give  some  indication  of  the  deficiencies  in  the 
measurement  and  data  reduction  assumptions.  Unfortunately,  many 
measurements  have  been  with  film  gages  which  have  doubtful  validity 
for  quantitative  turbulence  measurements,  since  substrate  thermal 
feedback  causes  probe  sensitivities  to  be  functions  of  frequency.  It  is 
particularly  serious  and  complex  for  multiple  films  mounted  on  the  same 
substrate  -  the  type  of  probes  used  for  shear  stress  and  normal  stress 

measurement.  Even  with  crossed-wire  probes,  data  interpretation  is 

involved  and  can  be  unreliable.  For  instance,  the  time-averaged 
expression  for  one  component  of  the  compressible  turbulent  shear  stress 
is  -(pv)’u'  whereas  the  hot-wire,  after  assumption,  measures  (pu)'v'  which 
differs  by  a  first-order  term.  Thus  it  is  clear  that  systematic 
investigations  of  fluctuating  velocities  are  still  needed  even  in 
zero-pressure  gradient,  compressible  boundary  layers  to  establish  a 
reliable  data  base  for  turbulence  modeling. 

The  effects  of  wind  tunnel  freestream  turbulent  flow  quality  are  known 
to  determine  model  performance  in  many  test  cases.  But,  for  decades, 
wind  tunnel  testing  has  mainly  been  conducted  in  test  section 

environments  which  have  not  been  adequately  or  consistently 

documented.  In  general,  the  effects  of  dynamic-flow  properties  on 
time-averaged  model  parameters  have  been  largely  ignored.  Perhaps  the 
major  and  most  widely  recognized  question  is  the  influence  of  freestream 


disturbances  on  model  boundary-layer  transition.  Recent  developments 
in  boundary-layer  transition  research,  particularly  those  of  the  NASA 
Transition  Study  Group,  have  stressed  the  dominant  role  that  freestream 
fluctuations  have  on  model  boundary-layer  stability  at  transonic  and 
supersonic  speeds.  Not  only  do  the  external  fluctuation  amplitudes 
influence  transition,  but  their  energy  spectra  are  also  particularly 
significant.  Streamwise  turbulence  also  produces  fluctuations  in  dynamic 
pressure  and  local  Mach  number  which  lead  to  time-dependent  inviscid 

forces  on  the  model.  The  normal  turbulence  components  produce 
fluctuations  in  the  angles  of  incidence  and  side  slip.  It  is  important, 
therefore,  that  we  document  the  dynamic  flow  quality  of  the  tunnels 

which  are  used  for  advanced  aerodynamic  testing.  In  this  way,  the  list  of 
tunnels  can  be  ranked  and  judgments  made  as  to  the  meaningful 
operating  ranges  of  adequate  flow  quality  in  each  facility  relative  to  each 

proposed  test  program.  Problems  of  wind  tunnel  flow  quality  have  been 

addressed  recently  in  ref.  3. 

The  purpose  of  this  report  is  to  address  the  requirements  and  pitfalls 
involved  in  the  use  of  hot-wire  anemometers  in  high-speed  flows  and  to 
detail  the  measurement  techniques  required  in  the  study  of  freestream 
and  shear-layer  turbulence. 
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Fig.  1.1  Reynolds  normal  stress  distribution  in  compressible 
turbulent  boundary  layers  (ref.  2). 


Klebanoff 


Fig.  1.2  Reynolds  shear  stress  distribution  in  compressible 
turbulent  boundary  layers  (ref.  2). 


Section  2  Definitions  and  Descriptive  Properties 
of  Turbulence 


Turbulence  is  the  usual  state  of  fluid  motion  except  at  low  Reynolds 
numbers.  It  is  a  three-dimensional  irregular  flow  condition  in  which  the 
various  properties  of  the  fluid  show  random  variations  in  space  and  time 
and  in  which  vortex  stretching  causes  the  fluctuations  to  spread  to  all 
wavelengths  between  a  minimum  determined  by  viscosity  and  a 
maximum  determined  by  the  boundary  conditions. 

Homogeneous  turbulence  is  statistically  independent  of  position  in  space 
but  is  not  necessarily  isotropic.  Isotropic  turbulence  is  statistically 
independent  of  direction,  ie.,  u'2  =  v'2  =  w'2  and  in  practice  is  also 
homogeneous  with  uV  =  v'w’  =  uV  =  0.  Scales  of  turbulence  are 
generally  referred  to  in  terms  of  the  dissipating  length  scale  (micro  scale) 
and  the  energy  containing  (macro)  scales.  Wave  number  is  the  inverse  ot 
wavelength  and  is  defined  as  the  frequency  divided  by  a  turbulence 
velocity,  normally  the  turbulence  convection  velocity. 

Mean  square  values  of  fluctuating  parameters  define  their  intensity. 
Probability  densities  provide  information  in  the  amplitude  domain. 
Information  on  turbulence  structure  in  the  space,  time  and  frequency 
domains  can  be  obtained  from  auto-  and  cross-correlation  and  power 
spectral  density  measurements. 

The  mean  square  value  of  a  random  property  is  defined  as 


L  M.A  akn  ut 


2.1 


■  litn  /Tu2(t)dt 
T-»«  1  o 

where  \|/x2  is  the  mean  square  value  and  the  square  root  \y  is  referred  to 

as  the  root  mean  square  (RMS)  value.  If  we  think  of  the  instantaneous 
value  of  a  fluctuating  quantity  as  a  time  mean  (steady)  and  a  fluctuating 

(dynamic)  component,  then  u  =  u  +  u'  where  the  mean  value  is  defined  as 

.T 

*  lim  /  u(t)dt  2.2 

T-»  °°  o 

and  the  variance  which  is  the  mean  square  value  about  the  mean  is 
expressed  as 

=  lim  /  Iu(t) -Mxl2cit  2.3 

T-*«  o 


that  we  have 


^ x  "  Mx  2.4 

where  the  positive  square  root  of  o2  is  called  standard  deviation. 

If  the  variable  of  interest  is  random,  it  can  only  be  defined  in  terms  of  its 
total  time  history  and  it  must  be  defined  in  probabilistic  terms.  There 
are  two  types  of  probability  density  function.  The  amplitude  probability 


function,  [P(x)],  defines  the  percentage  of  the  time  a  property  is  less  than 
a  given  threshold  level,  as  shown  on  the  next  page. 


It  is  given  by 

T 

Prob  [  x(t)  <  x  ]  =  P(x)  =  lim  -=£ 

T-*»  1 

The  complementary  distribution  function  is  defined  as  1  -  P(x).  The 
amplitude  probability  density  function  [p(x)]  defines  the  percen'.age  of 
the  time  a  property  has  a  value  within  a  given  window  width  as  shown 
below. 


We  define  the  probability  that  the  signal  will  be  in  a  given  window  as 
p(x)Ax  so  that 

p(x)  =  lim  Pr°b  1  *  <  x(t)  <  x*Ax  1  .  lim  lim  '  2  5 

AX-*0  AX  AX->0  T-*  <*>  1  AX 

ie.  p(x)  =  dP(x)/dx 


The  probability  density  distribution  of  a  continuous  variable,  which  is  the 
sum  of  a  large  number  of  independent  variables,  is  approximately 
normal  or  Gaussian. 

i©  p(x)  =  — — —  ©  2<f  2.6 

<fv2ir 

Most  random  processes  in  nature  are  Gaussian  or  have  related  Maxwell 
distributions.  Unfortunately  turbulence  is  not,  since  skewness  values  are 
not  necessarily  zero,  ie.,  u'2v'  *  0  for  instance.  Turbulence  then  is  the  sum 
of  a  large  number  of  processes,  but  they  are  not  quite  independent. 

Probability  densities  are  used  to  establish  a  probabilistic  description  for 
the  instantaneous  values  of  the  data.  They  are  shown  below  for  a  signal 
x(t). 


—  oo 


0 

Probability  distribution  of  x(t) 


+  oo 


Probability  density  of  x(t) 

The  auto-correlation,  or  co-variance,  describes  the  general  dependence  of 
values  of  a  property  between  one  time  and  another  and  is  defined  as 

1  * 

Rx(t)  =  lim  •=■  /  x(t)  x(t+t)  dt  2.7 

T-» »  *  o 

The  auto-correlation  is  always  a  real-valued,  even  function  with  a 
maximum  at  x  =  0  and  may  have  both  positive  and  negative  values  so 


Rx(-t)  =  Rx(t)  and  Rx(0)  £  1Rx(t)1  for  all  t 


Following  are  two  examples,  one  for  a  pure  sine  wave,  the  other  for 
turbulence  of  variable  frequency  content. 


Thus,  we  can  use  the  auto  correlation  to  separate  coherent  from  random 
data  and  get  an  idea  of  the  energy-frequency  content.  The  power 
spectral  density  function  Gx(f)  is  the  Fourier  transform  of  the 

auto-correlation  function.  It  can  also  be  determined  by  band-pass 
filtering  and  RMS  measurement  since 


2.9 


The  joint  probability  density  may  be  defined  as  before,  ie. 

p(x,y)  -  lim  lim  -&1 1 

Ax-»0  7  ^  T-»“  1  J 

Ay-»0 

when  x  and  y  are  statistically  independent. 

p(x,y)  =  p(x)  p(y)  2.10 

Cross-correlation  measurements  can  be  made  with  or  without  time  delay 

Rxv(t)  -  lim  4-  /Tx(t>  y(t+t)  dt  2.11 

y  T-»»  T  o 

It  is  not  necessarily  an  even  function  and  R  (t)  max  is  not  necessarily  at 

x  =  0.  But  we  can  say  R_„(-t)  =  R  (x)  Also 

*y  y* 

Rx(0)Ry(0)  2  IRxy(T)l2 
1/2[RX(0)  +  Ry(0)]  2  IRly(t)l 

If  Rxv  =  0,  then  x(t)  and  y(t)  are  uncorrelated.  If  R  =  0  for  all  times,  x(t) 
and  y(t)  are  statistically  independent  or,  if  the  time  average  *  0  then  R 

xy 

=  Wy 

Examples  are  shown  on  the  next  page  for  frozen  and  decaying  turbulence 
convection. 


no  decay 

n  n, 

time  delay 

As  we  will  see  later,  the  cross-correlation  is  useful  for  the  determination 
of  convection  velocities  and  determining  signal-to-noise  ratios. 

The  cross-spectral  density  (Gxy(f)]  is  the  Fourier  transformation  of  cross¬ 
correlation.  But  since  Rxy  is  not  an  even  function,  Gxy(f)  is  complex  and 
may  be  written  as 

Gxy(f)  =  Cxy(f)  -j  Qxy(f)  2.12 

ie. 

Gxy(f)  =  IGxy(f)l  e'j  0xy(f)  2.13 

where 

|Gxyl  -  /cx2y  *  QXy'  2.14 

6xy'  2»5 

'■'xy 

Typical  variations  of  these  parameters  with  frequency  are  shown  in  the 
next  figure.  A  useful  relationship  is 

GxGy  £  IGxyl^ 


We  generally  define  the  coherence  function  as 


When  Yxy2(f)  =  0  at  a  given  frequency,  x(t)  and  y(t)  are  incoherent,  ie. 
uncorrelated.  Both  x(t)  and  y(t)  are  statistically  independent  if  y*y2(0  =  0 
for  all  frequencies. 


Note  also  that  the  time  delay  x  =  0xy(f)/27if,  so  we  can  determine  the  time 
delay  as  a  function  of  frequency.  This  provides  the  relationship  between 
convection  velocity  and  the  turbulent  scale. 


Section  3  Turbulence  Measurement  Requirements 


Most  of  the  measurements  required  for  the  definition  of.  turbulent  flow 
fields  can  be  derived  by  considering  the  influence  of  turbulent 

fluctuations  on  the  overall  momentum  and  energy  balances.  Consider 
first  the  incompressible  form  of  the  streamwise  momentum  equation 
which  may  be  written  as 

3u  3u2  3uv  3uw  1  /  02u  32u 

at  ax  ay  az  '  9  ax  ^\ax2  ay2  az2/  3 

Now  if  we  write  u  =  u  +  u  and  noting  that  »  0 

andthat  -(G  *  u  )  *  -^-(u2  +  u72+  2  uu)  ■  2u-P-  +  |^- 

ax  3x  ax  ax 


we  obtain 


ax  ay  3z 


p 


ip 

3x 


+  ir 


(b2v  ,  a2 n 

\ax2  3y2 

du'w 


(  9u'2 

\3x2 


3uv 

_  JL.  - -  ■—  -f 

ax2  ay2 


az2 


a2n\ 

az2/ 

3-2 


The  last  three  additional  terms  have  the  units  of  stress  and  are  commonly 
known  as  the  turbulent  or  Reynolds  shear  stresses.  These  additional 
stresses  are  extremely  important.  If,  for  example,  we  have  a  pipe  flow 
where  the  velocity  fluctuations  are  ±10%  of  the  mean  flow,  then  the 
Reynolds  stresses  are  of  order  .OOlpu2.  Since  the  mean  velocity  gradient 
in  the  pipe  is  of  order  u/d  then  the  viscous  stress  is  |i(u/d).  The  ratio  of 
turbulent  to  viscous  stress  is  therefore  .001ud/p.  So  that,  for  a  pipe  flow 
with  a  Reynolds  number  of  100,000,  the  turbulent  stress  exceeds  the 
viscous  stress  by  two  orders  of  magnitude. 
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Thus,  it  can  be  seen  that  turbulent  shear  stresses  are  extremely 
important.  The  mathematical  difficulty  in  solving  the  equations  of  motion 
for  turbulent  flow  is  that  they  are  non-linear;  ie.,.  the  turbulent 
fluctuations  cannot  be  superimposed  on  the  mean  flow  without  affecting 
it.  Equations  for  first  order  terms  involve  double  correlations  and  so  each 
equation  involves  higher  order  terms.  The  problem  is  further 
complicated  by  the  fact  that  turbulent  cross-correlations  (stresses) 
depend  on  the  phase  as  well  as  the  magnitude  of  the  two  fluctuations. 


To  determine  the  turbulent  energy  balance  in  a  flow  we  need  a 
relationship  for  the  time  rate  of  change  of  the  individual  normal  stresses 


Djq2  Dyu'2  Djv’2  D^W2 
Dt  “  Dt  *  Dt  *  Dt 


where 


+  u  +  v  +  w  ~ 

dt  dx  By  B  z 


To  determine  this  relationship  we  substitute  u  =  u  +  u’  in  the  equation  of 
motion  and  multiply  each  term  by  u’. 

I  — £ 

Thus  ia .  u'is  .  U'|a' .  £I£l 

dt  dt  dt  dt 


u-n|a'.  U'2fa .  U'2|a'.u-n|a 

dx  dx  dx  dx 

U'2|a' 

dx  dx  dx 


Continuing  in  this  manner  and  collecting  terms  we  obtain 


I 

£ 

\4 


■  -  t  u'!£ *  lruVu' -(U'2|a .  uv|a ♦  u-vfs) 

'(U'*Jx  *  U  V  ly  *  u  w  |z')  » 

where  the  first  three  groups  on  the  RHS  represent  the  spatial  transport  of 
turbulent  energy  by  local  pressure  gradient  fluctuations,  the  viscous 
dissipation  of  turbulent  kinetic  energy  and  production  terms  ie.  rate  of 
mean  work  against  turbulence.  The  last  group  may  be  rewritten  using 


continuity  as 


-K 


3u,2v’  .  3u,2w' 


and  so  represents  the  spatial  transport  of  u’2/2  by  turbulence.  Thus,  we 
have  a  balance  for  the  turbulence  energy  in  a  flow  which  involves 
production,  dissipation,  and  convection,  where  the  production  term 
represents  an  extraction  of  energy  from  the  mean  flow. 

In  boundary  layer  flows  we  can  make  the  approximations 

3  _  n  3u  „  3u  .  n  - 
3  2  3x  3y 

so  that  the  continuity  and  Navier  Stokes  equations  reduce  to  the  following 


form: 


iS  ♦  ii  -o 

3x  3y 


niS „ vin. 
3x  3y 


J_  _3p  3  u  _  3u'V'  _  3u'2  ,  7 

P  3x  3y2  3y  3x 


Since  the  normal  mean  momentum  equation  may  be  written  as 


P  -  P  - 


where  p  +  — —  =  const 
•  2 


we  obtain 


a#  ♦via- 

3x  3  y 


u„4?«  + - y-(u'2  -  v'2)  3  9 

3x  dyi  3y  3x 


Thus,  we  can  solve  for  u  and  v  only  if  the  turbulence  terms  are  known 
and/or  can  be  modeled.  Normally  we  may  assume  that  (3/dx)(u’2  -  v'2)  is 
negligible  except  near  separation. 

For  turbulence  modeling  purposes,  it  is  usual  to  divide  the  boundary 
layer  into  two  regions,  an  outer  and  inner  layer.  The  outer  layer  contains 
80%  of  the  flow,  the  inner  layer  also  contains  the  viscous  sub-layer. 
Mixing  length  models  relate  the  local  shear  stress  to  a  mean  velocity 
gradient  by  means  of  a  turbulence  length  scale  or  mixing  length 


— r-:  T  2  I  3u 

uv  *  L 

1  3y 


310 


where  for  the  inner  region  L  =  Kjy  with  Kj  =0.4 

However,  in  the  viscous  sub-layer  close  to  the  wall,  we  must  account  for 
turbulence  damping  so  a  "damping  factor"  (A)  is  included  such  that 


L  .  K,y[l-  <*p(-£)] 


where  (A)  is  a  strong  function  of  local  wall  boundary  conditions 


ie. 
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a-2Mt)3 

so  that  (A)  approaches  zero  outside  the  sub-layer. 

In  the  outer  region  it  is  usual  to  assume 

-uV-K2f>0.4*||^|  3  13 

In  locations  where  boundary  layer  edge  intermittency  is  present  we 
multiply  this  expression  by  an  intermittency  factor 

1  '  erff ^ ("s’)'  0  7&1 

Y  =  - L-LLZ - L  3. 14 

2 

The  major  objection  to  the  mixing  length  hypothesis  is  that  it  is  based  on 
the  local  mean  velocity  gradient.  But,  we  know  that  turbulent  length 

scales  and  lifetimes  are  many  boundary  layer  thicknesses,  so  we  cannot 
assume  that  turbulence  properties  are  uniquely  related  to  local  mean 
profiles.  We  must  account  for  mean  flow  history  and  turbulence 

convection. 

Other  models  relate  the  turbulent  shear  stress  to  turbulent  kinetic  energy 
ie.  -u’v’  =  kq2,  where  k  is  taken  to  be  0.3.  This  assumed  proportionality 
once  again  has  only  limited  validity.  Along  the  axis  of  a  pipe,  for 

example,  the  stress  is  zero  but  the  turbulence  energy  is  not.  In  practice, 
therefore,  only  external  boundary-layer  flows  can  be  predicted.  It  cannot 
be  applied  to  wakes,  wall  jets  or  natural-convection  boundary  layers 
since,  in  these  flows  the  shear  stress  changes  sign  whereas  the  turbulence 
energy  cannot. 
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Thus  we  need  to  determine  the  types  of  measurements  which  will  be 
required  to  help  in  turbulence  modeling.  As  a  first  step,  let  us  determine 
the  mean  flow  kinetic  energy  equation  in  boundary  layer  form.  To 
achieve  this  we  must  use  u  =  u  +  u'  and  multiply  by  u.  So  that  for  the 


case  of  zero  pressure  gradient 

0ia  ♦  via  . 

3x  3y  3y2 


3-15 


becomes 


_  aiu: 

u — - — 
3x 


_ 

+  V  — — 

3y 


-t— i  3u 
u  v 

3y 


— -  ( u'v’u)  =  0 

3y 


3-16 


The  first  two  terms  represent  the  gain  due  to  advection  ie.  mean  energy 
transport.  The  third  and  fourth  terms  represent  the  loss  due  to 
turbulence  production  and  the  gain  due  to  energy  flux.  In  boundary 
layer  flow  there  is  a  loss  of  mean  flow  kinetic  energy  over  most  of  the 
layer  except  near  the  wall.  The  loss  in  the  outer  region  is  due  to  transfer 
to  the  inner  layer  by  means  of  the  Reynolds  stress  gradient.  For  flows 
with  pressure  gradient  the  term  u3(u'2/2)/3x,  a  term  dropped  in  the 
above  derivation,  cannot  be  neglected.  So  we  can  see  that,  in  general,  the 
mean  flow  energy  balance  involves  expressions  for  the  rate  of  change  of 
turbulence  kinetic  energy  and  shear  stress.  Derivation  of  these 
expressions  will  give  an  indication  of  the  required  measurements. 


To  derive  the  equation  for  the  rate  of  change  of  turbulence  kinetic  energy 
we  insert  u  =  u  +  u'  in  the  boundary  layer  momentum  equation,  multiply 
the  result  by  u'  and  take  the  time  average.  This  results  in 


3.17 


Dyq2 


1  /_  3q2  _  3q2 \ 


u'v'-|r-(l7lv'(p''‘!52)1*c) 


3y  v3y 


where  we  can  see  that  the  convection  of  turbulence  energy  by  the  mean 
flow  is  the  balance  between  the  turbulence  energy  produced  by  the  mean 


motion  working  against  the  turbulence  shear  stress  and  the  turbulence 


energy  diffusion  and  viscous  dissipation  (e). 


Finally,  the  rate  of  change  of  Reynolds  shear  stress  may  be  derived  by 
inserting  u  =  u  +  u'  into  the  boundary  layer  momentum  equation,  and 
multiplying  the  result  by  v'  and  taking  the  time  average.  Thus  we  obtain: 


Du'V 

Dt 


.2  3u  3  .  .2 1 

V  —  +  —  u  V 

dv  dv  J 


1  f  3p'u'  3p'v' 


3y  3y 


P  1  3y  3x 
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V2v’  +  V 


V2U'  ] 


where  the  rate  of  change  of  Reynolds  shear  stress  is  the  balance  between 
shear  stress  diffusion,  pressure  diffusion,  pressure  scrambling  and 
viscous  dissipation. 


If  heat  transfer  is  present,  by  analogy  with  the  mean  enthalpy  equation. 


we  determine  that 


Dje-2 


e'v2e 


J_  r  3u'6'2  3v'9'2 

2  L  3x  3y 


dve*  i 
3v  J 


1 

1 
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where  we  can  see  that  the  rate  of  change  of  enthalpy  function  depends  on 
the  conductive  diffusion,  production  and  turbulent  transport  respectively. 


By  analogy  with  the  Reynolds  shear  stress  equation  we  can  derive 

i 


Dv'e' 

Dt 


iE' 

ay 


+ 


v  -  - 

V  V2 ©'  +  IT  0'V2v' 
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In  this  case  we  have  two  sets  of  molecular  dissipative  terms  and  two  sets 
of  production  terms. 


Eddy  conductivity  may  be  expressed  as 

kT  =  -  pcp  vW 

Be 

By 


so  that  the  turbulent 


Pr 


Prandtl  number  becomes 

_  39 

_  e  cp  _  u'v  By_ 
r  kT  v'8'  3u 
3y 


3  2  1 
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In  most  cases  the  turbulent  Prandtl  number  is  assumed  equal  to  unity 
inferring  that  heat  and  momentum  are  transferred  by  the  same 
processes.  However,  measurements  (ref.  4)  indicate  significant  variations 
can  occur  across  shear  layers.  A  curve  fit  to  these  data  suggests  the  form 
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I 

I 

PrT  =  0.95  [l  "  0  5  (j)]  3  23 

Now,  all  these  previous  equations  express  the  rate  of  change  not  the 
magnitude  of  the  turbulence  kinetic  energy,  shear  stress  and  turbulent 
heat  transfer.  The  local  magnitude  depends  on  the  integral  along  the 
length  of  a  streamline.  To  get  some  idea  of  the  length  scales  of  turbulent 
motion,  we  can  study  spatial  and  space-time  correlations. 

From  the  turbulence  modeling  viewpoint,  information  on  the  turbulence 
scales  and  lifetimes  are  also  of  crucial  importance.  Since  turbulent  flows 
vary  not  only  in  time  but  also  in  space,  their  investigation  must  involve 
an  examination  of  both  the  spatial  and  temporal  statistical  structure. 

* 

>  Space-time  correlations  can  make  a  contribution  to  this  study  since  they 

give  evidence  of  the  heredity  and  structure  of  turbulence,  as  well  as 
values  of  the  convection  velocities  of  the  vorticity  and  entropy  modes 
compared  with  the  average  mass  transport  velocities.  Examples  of  both 
auto  and  space-time  correlations  in  a  compressible  turbulent  boundary 
layer  are  given  in  figs.  3.1  and  3.2.  These  data  were  obtained  on  a 

cone-ogive-cylinder  model  in  the  Ames  3.5-ft.  wind  tunnel  (ref.  5).  Fig. 
3.1  shows  the  auto-correlation  of  the  fluctuating  signals  on  the  cylindrical 
portion  of  the  model  176  cm  from  the  cone  apex,  at  two  positions  in  the 

turbulent  boundary  layer  and  in  the  far  field.  It  can  be  seen  that  there  is 

a  marked  variation  of  energy  distribution  with  frequency  across  the 
boundary  layer  and  that,  as  expected,  the  far  field  contains 

proportionately  much  less  energy  in  the  high  wave  number  range  than 
the  wall  region.  The  results  of  a  series  of  filtered  (4  kHz' 
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cross-correlation  measurements  at  several  wire  separation  distances  in 
the  boundary  layer  are  shown  in  fig.  3.2.  It  can  be  seen  that  each 
cross-correlation  curve  reaches  a  maximum  at  a  non-zero  value  of  the 
time  delay,  clearly  indicating  the  presence  of  convection.  The  amplitude 
of  this  maximum  is  a  function  of  the  wire  separation  distance.  A 
convection  velocity  of  these  disturbances  may  be  determined  from  the 
time  delay  at  which  the  maximum  of  a  particular  cross-correlation  occurs. 

The  peaks  of  the  cross-correlation  obtained  for  various  values  of  wire 
separation  distance  represent  the  auto-correlation  in  a  reference  frame 
moving  with  disturbances.  They  are,  therefore,  a  measure  of  the  lifetime 
of  the  disturbance  pattern  as  it  is  swept  along  with  the  mean  flow.  The 
long  turbulence  lifetimes  which  can  be  inferred  from  these  space-time 
correlation  measurements  (ref.  5)  illustrate  a  major  objection  to 
turbulence  models  based  on  local  mean  flow  gradients.  It  cannot  be 
assumed  that  turbulence  in  uniquely  related  to  local  conditions,  and  flow 
history  must  be  considered,  especially  when  attempting  to  calculate 
non-equilibrium  flows. 

To  conclude,  we  can  list  some  of  the  measurement  requirements  for  the 
definition  of  turbulent  flow  fields.  These  requirements  include: 

1)  Measurements  of  the  spatial  and  temporal  distributions  of  turbulence, 
kinetic  energy  and  shear  stress. 

2)  Determination  of  the  rates  at  which  these  properties  are  produced, 
transported  and  dissipated. 

3)  Determination,  from  spectra  or  correlation  measurements,  of  the 


contribution  of  different  length  scales  to  the  turbulence  kinetic  energy 
and  shear  stresses. 

4)  The  rate  at  which  Reynolds  stress  and  turbulence  kinetic  energy  are 
transferred  from  one  range  of  eddy  sizes  to  another. 

5)  Components  of  the  viscous  dissipation  from  point  time  histories, 
assuming  Taylor's  hypothesis  holds,  ie. 


TIME  DELAY,  msec 

Fig.  3. 1  Auto-correlations  across  a  compressible 

turbulent  bundary  layer  (ref.  5). 
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3.2  Space-time  cross-correlation  measurements 

in  a  compressible  turbulent  boundary  layer  (ref.  5). 


Section  4  The  Measurement  of  Turbulent  Fluctuations 

A  hot-wire  can  be  assumed  to  be  a  cylinder  placed  in  a  flow  field  for  the 
purpose  of  making  certain  types  of  flow  field  measurements.  The  wire  is 
assumed  to  have  length  L  and  a  diameter  d.  Furthermore,  we  assume 
that  the  wire  temperature  is  Tw  and  if  we  take  a  heat  balance  we  note 
that 

W  -  Qr-  Qc  -  Hw=  0  4.1 

where  W  =  the  energy  (electrical)  put  into  the  wire. 

Qr  =  Heat  loss  due  to  radiation. 

Qc  =  Heat  loss  due  to  conduction. 

Hw=  Heat  loss  due  to  convection. 

If  we  assume  that  there  is  no  heat  loss  due  to  radiation,  or  conduction  to 
the  supports,  then  equation  4.1  can  be  expressed  as  W  =  Hw.  In 
incompressible  flow  numerous  hot-wire  heat  transfer  studies  show  that 
Nu  =  Xd/kf  =  0.42Pr0,2  +  0.57Pr033Re°-5  4.2 

Equation  (4.2)  is  applicable  in  air  and  diatomic  gases  when  0.1  <  Re  < 
1000,  where  Re  =  p  ud/p  and  Pr,  p,  p  are  evaluated  at  the  film 

o  © 

temperature  Tp  =  (Tw  +  Tg)/2.  Free  convection  effects  may  be  neglected 

if  Gr  x  Pr  <  10‘4  and  Re  >  0.5.  For  a  5-micron  wire,  in  high-speed  airflows, 
Gr  x  Pr  =  106  and  the  Reynolds  number  based  on  wire  diameter  is 
generally  »  0.5. 

<■) 

Since,  the  heat  transferred  from  a  wire  of  length  L  is  X7idL(Tw  -  Tg)  =  I~Rw, 
we  may  rewrite  equation  4.2  as 


jikjL(Tw  -  Tg)[0.42Pr0-2  +  0.57Pr°-33Re°-5]  =  I2Rw 


4.3 


Now,  Rw  =  Ro[l  +  a(Tw  -  To)  +  y(Tw  -  To)2  +  ...  ]  and  neglecting  higher 
order  terms,  Tw  -  Tg  =  (Rw  -  Rg)/aRo  so  that 

I2Rv  -  [0.42Pr°2+  0.57Pr0  33Re°  5]  4.4 


ie 


where 


and 


I2Rv 
Rw  -  Rg 


A  +  B  /u 
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A 


0.42 


■rrk^L 
oc.  Ro 


Pr 


0.2 


4.6 


B 


0.57 


irkjL 
a  Ro 


4.7 


We  can  determine  A  and  B  experimentally  by  plotting  I2  vs.  'fa  for 
constant  electrical  resistance  as  shown  below. 

,t 


Now  let  us  consider  the  unsteady  heat  transfer  from  a  wire  in  a  turbulent 
flow  field  where  we  have  u  =  u  +  u',  Rw  =  Rw  +  rw.  If  the  sensor  is  made 
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N 


part  of  the 

constant 

current  circuit 

system  shown  on  the 

next 

page. 

assume  u'/u 

«  1  so 

that  rw/(Rw 

-  Rg)  «  1 

equation 

4.5 

may  i 

written  as 

I2Rw 

Wi.  rw  ) 

=  A  +  B  /ff  | 

4.6 

Rw  -  Rg 

\  Rw, 

)\  Rw-Rgy 

V  2uj 

Constant  Current  System 
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ire 


Rc  »  Rwire  SO  that  I  = 


Rc  +  Rwire  Rc 


or  since 


I2Rw 
Rw  -  Rg 


A  +  B  /u 


4.9 


we  have 


I2Rw 
(Rw-  Rg) 


,  *  B  /U  — 
2  2u 


so  that 


©  =  I  rw  *  - 


(Rw-  Rg) 


- b  ya-jt 

2IR,  U 


4.10 


le 


v/e^  = 


where 


28 


•  «*- 


and  e  =  -su.  Equation  4.10  applies  only  in  absence  of  thermal  inertia  of 
the  wire,  cooling  effects  of  the  supports  and  non-uniform  velocity 
distribution  along  the  wire. 

Let  us  first  consider  the  effects  of  thermal  inertia  and  rewrite  the  energy 
balance  to  account  for  the  heat  capacity  of  the  entire  wire,  C  ,  as 

I2Rw  =  (Rw-Ra)(A  +  B/u)  +  f>wcw  4.11 

where  cw  is  the  specific  heat  of  the  wire. 

Perturbing  this  equation  and  cancelling  the  steady  state  terms  gives 
I2rw  =  (A  +  B/u)rw  +  (Rw-R*)B>/U§  ♦  IT  412 

where 

c.  4.13 

or 

where 

Cw _ 

cc*R*[(A  +  B/U  )-I2 


M 


4.15 


4.16 


0(t)  = 


oc  gRg(Rw-  Rg)B/u 


2  C 


w 


But  from  the  mean  heat  balance 

A  +  B/U-I2  =  -R*  I2 
Rw  ~  R  j 


therefore 


Cw  (Rw  -  R g) 

I2  a  gRg 


4.17 


4.16 


The  solution  of  equation  4.14  is  obtained  as  follows.  Let  (prw)  =  X(t)  and 
prw  =  J  l.«  X.(t)dt 

-  6  liir  r"f  - x(t>  419 

-  7$*s  x(t) ’ *(m  lose’b 

.-.  B  - 

t  .(±)  JS) 

a  rw  *=  /  e  vn%(t)  e  w  dt  4.20 

-  go 

=  exp[-^  (t-t)]  dt 

Thus  if  u(t)  =  u*elon  where  u*  represents  a  velocity  amplitude,  then 

<Kt)  =  $*ei(a5t'n)  4.21 

where 

„  aRo(Rw-Ra)B /0  u*  (Rw- R*)2  B /U  u* 

**  = -  =  - - -  4.22 

Gw  2u  MI2Rft  2u 
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Equation  4.20  then  gives 


rw 


.  iut 


**  M 


1  +  iwM 


rv,.ei<o«-y> 


4.23 


where 


rw 


-  g*  M 
>/l  +  w2M2 


and 


^  -  tan'V+wM) 


If  we  take  e  =  -su  where  s  is  the  velocity  sensitivity  then  we  can  see  that 
thermal  inertia  affects  both  the  amplitude  and  phase  of  the  hot-wire 
response  ie. 


s  » 


J 1  +  <^2M2 


‘NTI 
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where  NTI  denotes  no  thermal  inertia.  This  shows  that  both  phase  shifts 
and  amplitude  variations  are  functions  of  the  frequency  of  the  velocity 

fluctuation.  M  is  defined  as  the  time  constant  and  coQ  =  1/M  is  called  the 
roll-off  frequency.  The  time  constant  is  a  strong  function  of  wire 
diameter  and  mass  flux.  Take,  for  example,  a  5  micron  tungsten  wire  in  a 

turbulent  flow  of  30  m/sec,  with  Tw  -  Tg  =  115°C,  Tg  =  288°K,  and  I  =  76 
mA,  the  wire  time  constant  is  approximately  3  x  10‘4  sec  and  the  signal 
amplitude  at  250  Hz  is 

s  _  1  =  1  _  1 

SNTl  J \  Z  9  x  10'8(2ir*250)2  / 1  +  2  2  1.1055 

so  that  the  signal  has  already  dropped  down  10%  at  very  low  frequencies. 
The  signal  phase  angle,  which  may  be  calculated  from 


<(>  =  tan'^-coM)  =  tan'^rc  x  250  x  3xl0'4)  =  tan_1(0.47)  4.25 

shows  that  there  is  a  25-deg.  phase  lag  at  250  Hz.  Thus  it  is  crucial  that 
we  compensate  for  these  effects. 


! 


A  typical  compensation  circuit  is  shown  below 

C  Re< 1  = 


RiZc 

Ri+Zc 


.1 


R2 ;;  vo 

jj 


An  analysis  of  this  circuit  gives 


R*q+  R2 


Vo  = 


V1R2 

Req+R2 


=  (i  *  — — ^ 
K2  V  RlR2  +  R2Zc/ 


4.26 


4.27 


Now  again  assume  u  =  u  +  u'  ,  Rw  =  Rw  +  rw  ,  Tw  =  Tw  +  tw  and  since 

Rw  -  R  4 

rri  w  /  a  a 


Tw  -Tg  = 


Rft 


4.26 


so  that  rw  =  Roatw.  On  substitution  in  equation  4.11  we  obtain 

1 2 Rw  + 1 2 rw  =  (Rw+ rw~Rg)[ A  +  Bv/iJ  ^  1  +  ^  =  4.29 

drw 


cci'Rg  dt 


Since 


I2Rw=  (a  +  B/u)(rw-R*) 


4.30 


neglecting  higher  order  terms  we  obtain 


VJVr1 


V 


*.*l1 


R2  /  1  +  oo  Rl  C 

Rl/  (.  R2\2  2d2  2 

y  ii+ rt)+ 00  Rzc 


j  tan1  r 


4.36 


Clearly  then  we  can  match  for  the  hot-wire  amplitude  roll  off  if, 


1  +oo2M2 


/  1  ♦  oo2  Rl  C2 

(l  +  |i^+  r!  c2 


4.37 


Now  defining  to'  =  1/RiC  and  to  =  1/M  for  the  hot-wire,  for  amplitude 


compensation  we  have 


2  1+(^) 

i.(*-)  - - 

/  D-»  \“  /'D'> 


4.36 


and  the  circuit  phase  angle 
/  R2  00  \ 


oo_  _  /  Rl  ^0 

V  1  R2 

V  1  +  TTT 


/ 


4.39 


for  hot-wire  phase  correction. 

In  an  actual  compensating  circuit,  the  capacitance  must  be  a  variable  as 
the  wire  time  constant  depends  on  both  overheat  and  mass  flow  (pu). 
Usually  R2  «  Rj  (say  1/1000)  so  that 


mwmsimmsm 


**  n  Vo  ~ 

at  "  ’  0  Vf  *  Rl 


4.40 


and  for  w  -+  »  — — >  1 

Vi 


The  point  on  the  curve  where  Vo/Vi  flattens  off  is  cd/cd  =R2/Ri  or  to 


=  l/CR2so  that  variable  capacitance  provides  for  adjustable  limiting 
frequency  compensation  but  the  overall  gain  is  fixed  at  R2/Ri.  For 


variable  R  l,  the  frequency  at  which  compensation  is  correct  is  fixed  at 


comp 


=  1/CR2  and  the  overall  gain  changes  with  Ri.  Since  it  can  be  shown 


y'  is  maximum  when 


That  is  when 


e)  - 


.  R2  R2  / irj  \2  ~ 

1  +  rT  rI  (-J  =0 


—  =  7i  +  ^ 

<*a  y  R2 


i  V,-^> 


G  =  gain  factor  and  coA  =  set  for  MA 


Since  the  wire  output  may  be  written  as 


where  =  1/co 
o  o 


4.44 


4.45 


for  co  ->  0  Vo/Vf->  G 

(o  — >  °°  Vo/Vf  — »  G  o>o/co A 

We  can  see  from  figure  4.1,  that  as  long  as  coo/coA  is  within  ±  30%  of  unity, 
the  response  is  flat  above  gj/g>0  *  5.  So,  apart  from  signals  of  frequencies 
less  than  5coQ,  we  can  account  for  errors  due  to  incorrect  circuit 
time-constant  settings  (coA)  by  multiplying  the  signal  by  coA/coo. 


Alternatively,  we  can  write 

ef  -  X~F(f)  df  4.46 

where  eo  is  the  uncompensated  wire  output  and  F(f)  is  the  power 

spectrum  ie.  the  signal  amplitude  squared  vs.  frequency.  The 
compensated  signal  is  then 

eJomp  */0*F(f)(l  +  w2+  Mi)  df  =  e2  +  2F(f)  df  4.47 


The  true  signal  would  therefore  be 

-/*F(f)(l +oo2+MT)df  =  e?  +  Mt  f  w2F(f)df  4.46 
1  0  0 

eliminating 

/  ^2p(f)  df 
0 

we  have 

ecomp~  gp  Mt  ^ 

e|-e*  Me 

which  may  be  rewritten  as 

eT  =  Ml  L  ^|-ec2omt>) 

ecomP  Me  yj  e|omp(e|-eJ) 


4.49 


4.50 


This  correction  is  acceptable  when  Mt  and  Me  differ  by  less  than  30 
percent. 
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Fig.  4.1  Hot-wire  compensating  circuit  response. 
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Now  let  us  turn  our  attention  to  the  constant  temperature  method  (rw  -» 
0)  and  consider  a  hot-wire  probe  in  the  circuit  shown  below. 


Now  assuming  1  =  1  +  i,  u  =  u  +  u’  and  substituting  in  equation  4.5, 
neglecting  terms  in  rw,  we  obtain 

I2R  +  2iIRw  -  (Sw-R4)(A  ♦  B/tT)  +  (Rw-R*)B/u  ^  4.5 1 


so  that 


— v:Ra  B/ir 

4lRw  G 


4.52 


and  e  =  iRw  =  sct  u.  Where  the  constant  temperature  sensitivity  is  given 
by 

(Rw  -  R$)B/u 


>ct 


4l0 


453 


We  also  see  that  the  constant  temperature  sensitivity  is  related  to  the 
constant  current  sensitivity  (s)  by 
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so  that  for  Rw/Rg  >  3/2,  s  <  s,  ie.  at  low  overheats,  constant  temperature 
operation  is  less  sensitive  than  constant  current. 


In  order  to  study  dynamic  behavior,  we  must  consider  small  changes 
(lags)  in  Rw  =  Rw  +  rw.  Neglecting  higher  order  terms  in  i,  rw  and  u',  we  see 
from  the  thermal  equilibrium  equation  that 

11  JT 

T2rw  +  2ilRw  =  (A  +  B/u)rw  +  (Rw  -Rg)  B/u  r=  +  Gw— w  4.55 

2u  at 


using  Tw  =  rw/aRo  and  substituting  for  gtr  we  obtain 

di  ocRo  ,  _ ,  _  o  _  n —  occy+rl  Ro(Rw  -  R  g)  ti 

+  (A  +  B/1T  -  1 2  +  21 2Rw  g,r )  i  =  ■  B'r  „ - B/u 


dt 


Gw 


2u 
4.56 


solving  this  equation  as  previously,  we  determine  that 


Met 


_ Gw _ 

ocRo  (A  +  Bx/vT  -  I2  +  2l2Rw  gtr  ^ 


4.57 


or 


Met 


M 

1  +  2 

Rw  -  R  g 

Rwgtr 

R* 

4.56 


ie.  Met  «  M. 

Typically,  Rw/Rg  =  2,  Rw  =  10  ft  and  g(r  =  10  mhos,  so  that  Met  =  1/200  M. 
However,  we  also  note  that  the  response  depends  on  (Rw  -  Rg)/Rg,  so  that 
when  Rw  -  Rg  0  Met  M.  Thus,  the  constant  temperature 
anemometer  frequency  response  is  poor  at  low  overheat  (ref.  6). 


The  phase  shift  is  also  a  strong  function  of  overheat  and  mass  flow  rate. 
At  low  overheat,  phase  shifts  in  the  feed  back  amplifier  can  lead  to  large 
bridge  oscillations.  To  avoid  this  we  need  flat  amplifier  response  to  at 
least  twice  the  desired  compensation  frequency.  Bridge  stability  also 
depends  on  the  amount  of  out  of  balance  which  is  a  function  of  Rw  -  Rg. 


Typically  for  a  5pm  tungsten  wire  at  high  overheat,  -3db  relative 
amplitude  can  exceed  100  kHz,  but  with  a  phase  lag  of  60  deg. 


Now  let  us  consider  the  effects  of  conduction  to  the  wire  supports. 
Generally,  hot-wire  supports  are  much  thicker  than  the  wire  itself  for 
reasons  of  strength  so  that,  since  their  relative  resistance  is  low,  they  will 
not  be  heated  appreciably  by  the  electric  current.  We  can  determine  the 
hot-wire  temperature  distribution  and  determine  the  effects  of  heat 
conduction  to  the  supports  by  rewriting  the  energy  balance  as 


I2Rw  =  (Rw  -  Rg)(A  +  B/u)  -  —  kw 


d2Tw 


4.59 


where  kw  is  the  thermal  conductivity  of  the  wire.  Now  since  Rw  -  Rg  = 


aRo  (Tw  -  Tg)  we  obtain 

d2(Tw-Tg)  aRo[A  +  B/u-I2l 


I2Rg 

(Tw-Tj)t— 5 —  '0  4.60 

24-  kw 


where 


Tw  -Tg  *  0  for  x  =  i  ~2 

d(Tw-Tg)  ^  . 

-  =  0  at  x  =  0 
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The  solution  of  this  equation  is 


T w  -T*  = 


T2^i 

1  Ro 


oc  (A  +  B/u-I2) 


1  - 


cosh  (x/Lc) 


4.61 


cosh  (x/2Lc)J 

where  Lc,  the  cold  length  ie.  the  length  effectively  cooled  by  the  end 
supports,  is  a  function  of  wire  diameter  f(d)  and  may  be  written  as 

4.62 


L  -  — 

i-c  -  — 


irkw 


2  J  oc  (A  +  B/u  - 1 2 ) 

and  the  effective  wire  length  is  equal  to  1 

approximations  equation  4.61  may  be  written  as 

cosh  x/Lc 


Tw  -Tg 


=  1  - 


Tw^-Tft  “  cosh  L/  2LC 

and  examples  are  shown  below. 


2(Lc).  With  some 


4.63 


oo 


Fig.  4.2  Hot-wire  temperature  distributions. 


From  this  figure  we  can  see  that  we  must  keep  L/2(LC)  >  5  ie.  no  more 
than  about  20  percent  of  the  wire  cooled,  no  more  than  10  percent  at 


each  end.  For  typical  operating  conditions  this  means  L/d  =  100  to  200 
for  platinum-rhodium  wires  and  still  higher  for  tungsten  since 
kw(tungsten)  =  2.5  kw(platinum). 

The  effects  of  end  losses  on  dynamic  response  are  very  complex.  They 
may  be  summarized  as  follows.  It  is  impossible  to  correctly  compensate 
with  simple  RC  circuits  since  the  phase  shift  is  not  constant  along  the 

wire.  However  deviations  are  small  (within  2  percent  for  L/2(Lc)  >  3)  if 

compensation  is  correct  for  co  =  0  and  °° .  The  wire  time  constant 

decreases  as  L/d  increases,  but  there  is  a  limit  to  L/d  due  to  strength  and 

non-uniform  velocity  distribution  along  the  wire  which  we  will  consider 
next. 


Consider  the  case  of  uniform  flow  u  and  turbulence  intensity  up  Since 
the  wire  voltage 


e  =  Kf  u.(xjdx. 
o  1  2 


4.64 


then  if  u,  =  constant, 

e*-K2L2u?  4.65 

This  would  also  be  the  case  in  non-uniform  flow  if  we  measure  a  true 
"point"  value.  Now  if  R2(x2)  is  the  lateral  correlation  coefficient 

distribution  along  the  wire  then,  when  the  two  rms  levels  u'1(x2)  and 
u'j(x2+  a)  are  equal. 
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4.66 


U!(x2)  ut(x2+a)  _  u1(x2)  u^x^a) 
R2  *2  ’  ui(x2)ui(x2+a)  "  ui(x2)2 


But  if  u’j(x2)  is  not  constant,  the  measured  mean  square  voltage  output  is 


em=  K2[  jTI'ui(X2)dX2]4 


4.67 


which,  if  Jf(x-)  =  |f(x2)  can  be  written  as 


e^=K2jfjf  u,(x2)  u1(x'2)dx2dx2 


4.66 


In  homogeneous  turbulence 

u2  =  u2(x2)  -  u2(x2) 


so  we  have 


em=  r2  u?/oL/oLR2(x)(x2-x2)  dx2dx2 


4.69 


R2(x)u? 

0l<X'2)  -  T~7 - f 

Uj(x2) 


4.70 


Now,  substituting  in  equation  4.68,  and  defining  s  =  (x2  -  x2’)  we  obtain 


e2  =  2K2uf/L(L  -  s)R2ds 

m  jQ 


comparing  e2  and  e  2  given  by  both  equation  4.65  and  4.71  we  see  that 


e2  T  2 
-==■  =  tm 
e2 
m 


where 


/L(L  -  s)R2  ds 
Im  L2  0 
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Since  R2  <  1.0,  the  turbulent  energy  measured  with  a  long  wire  is  always 

smaller  than  the  true  point  value.  So,  to  correct  em2,  we  need  to  know  the 

lateral  cross-correlation  coefficient.  Now,  if  the  wire  length  far  exceeds 
the  lateral  integral  scale, 

A  ■  Jf  R2  dx2 )  4.74 

then 

2K2u|(-6A2-A2>  4.75 


depends  on  the  shape  of  the  correlation  curve.  When 

~ 0  so  that  -♦  0 
A2 


ie.  no  turbulence  would  measured  at  all  if  L  =  °°. 


The  effects  of  L,  d,  and  L/d  can  be  summarized  as  follows.  For  a  small 
time-  constant,  ie.  high-frequency  compensation,  we  require  a  small  wire 
diameter.  For  a  long  uniform  wire  temperature  and  high-frequency 


response,  we  require  L/d  >  200.  Decreased  wire  diameter  increases 
signal/noise  ratio,  but  decreases  strength  and  large  L/d  probes  give  poor 
response  in  a  variable  field.  However,  small  L/d  probes  lead  to  adverse 
aerodynamic  interference  effects.  So  we  can  see  that  probe  design  must 
be  a  compromise. 


The  effects  of  mean  fluid  temperature  changes  can  be  assessed  as  follows. 
Consider  the  energy  balance 


I2Rv 
Rw-  R g 


ACT*)  +  B/u 


4.77 


where  A  =  f(Tg)  whereas  B  is  generally  constant  provided  Tw  is  adjusted 
to  keep  Tw/Tg  constant.  But  the  turbulence  measurements  are 
independent  of  A  since 


I2  Rwfw 
(Rw-  Rg) 


4.76 


so  that  mean  temperature  changes  do  not  effect  hot-wire  turbulence 
measurements.  The  primary  effect  of  mean  flow  temperature  changes  on 
the  wire  calibrations  are  shown  in  fig.  4.3. 


Nu 


T* 

•zr  =  constant 
Tg 


Fig.  4.3  Effect  of  ambient  temperature  change  on  hot-wire  calibration 

Finally,  let  us  consider  the  effects  of  large  turbulence  fluctuations  and  the 
limitations  they  impose  on  our  earlier  assumptions.  In  such  cases  the 
effective  velocity  may  be  written  as 


u 


•11 


1  “i  .±?l\±h\  1  u*  3  “.“I  , 


2  a  I6jl  }  j! 


4.79 


If  we  time  average,  then  the  measured  velocity  is  given  by 


rz — '  r=>  HZ'  r  1  u-2  1  U-2  1  u,3  3  u.u2  -1 

^  [1_65"2  +  4  ia+T5is'6”5r+  ]  4  60 


r  i  u:  \  \iy  i  u/  3  u.u‘  i  ,  . , 

••  V,.  L1  4  u2  2  u2  6  u3  4  u3  J 

So  if  the  fluctuations  are  large,  the  measurement  is  a  function  of  u,,  u2 
and  Uj u 2  plus  higher  order  terms.  So  we  must  correct  the  mean 
measurements.  Neglecting  higher  order  terms  we  have 

r  l  u2  i  u  2  -I 

Jail 
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In  isotropic  flows  Uj2  *  u22  so  that  umeas  >  uactual.  Now  the  hot-wire 
response  to  turbulent  fluctuations  may  be  written  as 


I2Rw 


Rw-  Rg 


putting  Rw  =  Rw  +  rw  where  rw  is  not  «  Rw.  We  obtain 
I2  (Rw+rw)  =  (A  +  B\/ru')(Rw-Rj+rw)  +  B\/U'(Rw-Rj+rw)  u*  4.63 


where 


ill’ll  1  3  , 

u  ‘  2  n  8  u2  4  u2  16  u3  8  u3 


4.64 


averaging  we  obtain 

I2Rw  =  (A  +  B/u)(Rv-R6)  +  B/u(Rw'Rg)u*  +  B/u  u*rw  4.65 

Now  for  linearized  theory 

I2Rwlin=  (A  +  BnATMRw^  Rj)  4.66 


so  that  from  equations  4.83,  4.85  and  4.86  we  obtain 

/t>  n  \  2  O 


_  (Rw«i-Rj)  u2  r  l 

e2  =  I2  rw  =  (Rw-R*)2  — Tjr - B2  u  — 1—2  [l  -  0(Uj,u2)J  4.67 


"  1  1W  v . I  Rj  4u2L  2  J  v 

According  to  linearized  theory,  equation  4.10,  we  determine  that 

4.66 


or  after  further  calculation 

m  m 


(M;  OU1  •’‘(U‘-U2>1 


4.90 


where  \ji(u1,u2)  is  a  function  of  Ujn,  u2m,  Ujnu2m.  So  that,  in  highly 

turbulent  environments,  what  is  actually  measured  is  not  the  u; 

component  alone  but  a  mixture.  To  measure  Uj  we  must  not  only  know 

u2  but  also  the  higher  order  correlations  which  are  not  usually  known  for 

non  isotropic  turbulent  flows.  We  may  conclude  that  the  correct 
measurement  of  turbulence  of  high  intensity  is  practically  impossible 
with  a  hot-wire.  But,  if  errors  of  10-20%  are  accepted,  it  is  permissible  to 
use  linearized  theory  up  to  turbulence  intensities  of  20-25%.  This  covers 
a  large  range  of  flows  as  the  table  below  indicates.  Errors  due  to  high 
turbulence  and  resulting  directional  intermittency  are  discussed  in  ref.  7. 


Table  1 

Flow  Type 

Wind  Tunnel  (poor) 

Wind  Tunnel  (good) 

Turbulent  Jet 
Karman  Vortex  Street 
Turbulent  Boundary  Layer 
Turbulent  Wake 

Screen  Turbulence  (nearly  isentropic) 


Turbulence  Intensity 
0.1  -  0.5% 

<  0.05% 

30  -  60% 

5  -  15% 

4  -  10% 

2  -  5% 

0.3  -  2% 


The  measurement  of  turbulence  characteristics  with  a  hot-wire 
anemometer  may  be  achieved  as  follows.  Consider  the  response  of  the 
hot  wires  designated  0,  I  and  II  shown  in  the  flow  on  the  following  page. 


I  SI  \ 


I  /  Ui  sm  P 


V  =  Uj  sin  P  sin  0 


For  thermal  equilibrium  we  may  write 


eo  =  *  (si  )0U  i 


ei  =  -  <Sl)lUl  '  <S2>IU2 


4.92 


e!I  *  ^S1  ^IIU  1  +  ^S2^IIU2 


4.93 


Now,  if  we  use  an  X  array  with  identical  wires  or  make  the  sensitivities 
equal  by  adjusting  the  individual  currents,  we  get 


(Sj)j  —  (Sj)jj  and  ( ) j  —  (Sj)jj 


if  pj  =  p2  =  45  deg.  then  Sj  =  s2  =  s,  so  that  Uj,  u2  and  UjU2  may  be 


obtained  from 


«i-<V*n)2‘4s*ui  '4siul 


*i>-(Ven>2‘4s2u!  “4s2ut 


4  94 


and  ex2-en2  -4  8^2-4  8^2 


where 


'l " 


^  «  const,  V  -  const 


4.95 


S2  * 


3e 

3V  v 


4.96 


The  calibration  sensitivities  can  be  determined  as  follows.  Since 


3e  _  3E/3v 
3v  3*/  3* 


4.97 


and  velocity  component  resolution  shows  that 


3V  dQ 

[sin  *  cos  8  ^  ♦  cos  0  sin  0] 


4.96 


If  the  wire  is  at  45  deg.  to  the  flow  then 


3V  U 

37  -  Ul 


4.99 


1  3E 


3v  lu  ui  d*  I  Uj  x  const 


4.100 


so  calibration  is  relatively  straightforward. 


Alternatively,  UjU2  may  be  obtained  by  measuring  the  correlation  factor 


between  es  and  eD  ie. 


^AB  =  ^AB<trvf8«)  9aPb  =  ^A^BeS  eD  =  U1U2 


where  G.  and  GR  are  the  anemometer  gains.  So  that 


4.101 


U.U,= 

1  2  4s2 


To  determine  the  cross-correlation  coefficient  of  Uj  we  put  two  wires  at 
points  A  and  B,  so  we  have 

eA  =  <S1>A<U1>A  4*102 

and 

eB  =  ^Sl^B^Ul^B  4.103 

The  cross-correlation  of  Uj  can  then  be  calculated  from 


<VeB>  "  (VeB)4 


4eAeB 


2  (s1)a(s1)b(u1)a(u1)b 


4.104 


^eA+  eB^  +  ^eA_  eB^  2  (  eA  +  eB  ^  ^S1^A  ^UPa  +  ^S1^B^U1^B 


If  we  adjust  =  eBz  ie. 

cross-correlation  coefficient  reduces  to 


eA  eB  " 


=  eR2,  then  the 


2eAeB 

eAeB 


2eAeB  _  2eAeB 


where  *  denotes  RMS  values  ie. 

TO. 


ReAeB  ■ 


<U1>A<U1>B 


(ut)‘ 


4.105 


We  can  determine  the  Reynolds  shear  stress  coefficient 


U1U2 

V2 


from 


-  ©2 
i  en 


4s2ulU2  _ 
4s2uju'2  '  uiu*2 


4.106 
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So  we  do  not  need  to  calibrate  the  probes. 


To  determine  the  space-time  correlation  of  the  Reynolds  shear  stress 

(Ui  u2  )(Uj  u2  )  =  4107 

A  A  B  B  A  B 

we  can  use  two  hot-wire  X  probes  and  measure  the  correlation  function 
between 

*»*  ‘ (V W '  2SaU‘2*  .... 


4.106 


%  =  (%’  enB)  "  2SbUi2b 


(V  V2-  (v  v2 
(v  %)2+(v  V2 


4  eDAeDp  2SASBU12  U12 

- £-2-  =  — = - 4.109 

2(V  W  SAU12  +  V*12 


V »A*  V  V  V  A  12a  B  12b 

If  we  adjust  eD^=  eD2B  ie  e*  Ae^  -  eD2A  *  eD2  then 


RUi2u12  = 

A  B 


'VDb  ba  -db 


V^B 


V^B 


**A% 

ei>Aei>B 


4.110 


so  that 


A  B 


U1 2  ^1 2 

U12aU12_ 

U12aU12d 

A  B 

A  B 

A  B 

i 

i  w 

1<M^ 

1  3 

2.U1 2 

A 

B 

A  B 

4.1 1 1 


Triple-point  correlation  measurements  determine  the  kinematic  diffusion 
of  turbulent  energy.  Terms  such  as 


U2U2 


may  be  obtained  from  an  X  wire  when  Sj  =  s2  =  s  as  follows,  with 


,2  _ 


=  (ej  +  eu)2  =  4s2  Uj2  and 


B  =  ej  -  en  =  2sir 


then 


(A2+  B)2  -  (A2-  B)2 


4  A2B 


2  s3u5u2 

( A2+  B)2  +  ( A2-  B)2  2 ( A4+  B2)  s4u}  +  s2u| 


4.112 


If  we  adjust  A4  *  B2,  then  A2B'=A4=B2<  so  that 

ufu2 

R  =  — — — 
u'|tf2 


4.113 


or  we  can  measure  the  correlation  factor 
<eA+ 


so  that 

2  AB  -  GaGb  RAIW)  *  Rab  .  AB  =  <114 

Then  divide  by  say  s3u'3  to  remove  the  wire  sensitivity, 
components  eg. 

u3  &  u2u2 

may  be  obtained  in  a  similar  manner. 


The  space-time  triple  point  correlations  can  be  obtained  from  two  ] 
at  two  different  points  in  the  flow.  Defining  A  and  B  as 


=  *2 


A  =  e 


i 2  _  p2  •  q2 


X  S*\\ 


B  -  el  ■  eiB-  ®fl/  SIU1bU2b 


4.1 15 


then 


Other 


wires 
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'  (U*)A(U1)B  ^  (U2)A(Uj)A(Ui)B 


4.1 16 


R(xt)s.(^v 

(uj  )* 2  (UjV  *A  18  "A  *  A  *  B 

A  A 

Flatness  factors,  which  give  an  indication  of  the  extent  of  turbulent 
intermittency,  are  defined  as 

—  =  *  T« 

u'«  [ual2 

V4  ,  w4 

—  =  Tv  &  -  =  Tw 

.4  .4 

V  w* 

The  turbulence  intermittency  may  be  determined  from 

v.JL.&.AL 

Tu  Tv  Tw 

Longitudinal  microscale  may  be  defined,  using  the  auto-correlation  Rn,  as 

2 


2Uj 


1  ©)2  [£}  r“(x‘-oao) 


4.117 


X£=  0 


©)-u‘feR“(x,Hxi-o  4114 

for  steady  homogeneous  turbulence  only.  Since  for  steady  homogeneous 
turbulence 

Rn(XiAO,0)  -»  l  -  x2/*2  as  r->  0 
Now  invoking  Taylor's  hypothesis  x  =  uct 

\  =  VUc 

where  is  defined  by 
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g 

V? 


-s 

ft 

i 


Ru( 0,0,0,t)  -»  1  -  — -  as  t-»0 


ft 


I 


1 


*w 

$ 

I 

ft 


m 


?? 

c*» 

•ft 

<!> 

$ 

i$l 

!/*( 

!« 


Is 

$ 

ft 

1 


ie.  X,  = 


4.119 


Similarly,  we  may  measure  the  lateral  microscales  defined  as 


2  2U> 

S'  7T— rs 


4.120 


fey 


fey 


4.121 


We  can  also  measure  terms  such  as 


fey 


which  appear  in  the  equation  for  the  turbulent  dissipation  (e),  by  placing 
two  wires  a  distance  x2  apart  to  obtain 

[(u^a  -  <Ui>b12 

With  a  differentiation  circuit,  which  will  be  described  later,  we  can  also 


obtain 


(trf 


and  in  a  similar  manner  we  may  measure 


w 


mm*\  w*  »m  w*  tP  mmm% 


m  *  © 


When  the  electronic  noise  is  large  enough  to  influence  measurements 
made  with  a  hot-wire  anemometer,  the  noise  must  be  substracted  from 
the  anemometer  output  to  obtain  the  correct  net  result.  To  make  this 
correction,  the  hot-wire  equation  can  be  expressed  (ref.  8)  as 


(e‘  +  e^)  (u‘  un) 

GE  u 


4.122 


where  the  gain  is  G  and  the  noise  is  denoted  by  the  subscript  n.  The  level 
of  the  electronic  noise  can  be  approximately  obtained  by  covering  the 
hot-wire  probe  and  measuring  e'n  and  En.  Under  this  condition,  equation 
4.122  can  be  expressed  as 


*  S 


(f) 


4.123 


Forming  the  mean  square  of  equation  4.122  and  assuming  that  there  is  no 
correlation  between  the  electronic  noise  and  the  velocity  fluctuations  we 
have 

r  -  -  _ 

.2 


te'  *  en)2  S2 


G2E2 


4.124 


So  that,  the  velocity  fluctuation  in  the  flow  is 


(£) 


(e'  +  en)2 


S  2G2E2 


-m 


4.125 


Substituting  equation  4.123  into  4.125  results  in 

(e1  +  eft)2 


(f)  - 


n  \  2 


s  2G2E2 


-GH 


0  's2g„2e* 


4.126 
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w:* 

$$ 

V*fr 


Aj 


m 


,1 

«8 


;?w 


»<W 


A  value  is  now  required  for  the  quantity  u  2/u2  in  equation  4.126.  From 


the  calibration  of  the  hot-wire,  we  have 


log  E  =  A  log  u  +  log  B 


4.127 


which  can  be  written  as  follows  by  noting  that  A  =  S 


« -  (if 


4.126 


Then  using  equation  4.128  un/u  becomes 


la  =  fJL^s 

n  / 


4.129 


Substituting  equation  4.129  into  4.126  gives  the  final  result 


u‘\2  (e'  +  e‘)z  /E  As  el* 


(})■ 


S2G2E2  VE  )  S2Gn2E2 


4.130 


or  in  terms  of  the  rms  voltages  measured  with  an  rms  meter 


2 


fHT  -  fEn\s  e„ 

Vu7  S2G2E2  'E'S2Gn2E2 


4.131 


where  the  subscript  T  refers  to  the  total  rms  measurement. 


The  effect  of  noise  on  correlation  measurements  may  be  assessed  as 
follows.  Let  the  noise  be  a  and  b  on  two  anemometer  outputs  v}  and  v7, 

then  the  cross-correlation  function  may  be  written  as 
„  (v.  +  a)(v,  +  b)  v7T2 

T>_  *  “  —  i  C  s.  ion 


((Vj  +  a)2  )'  ((v2  +  b)2  Y  ((vf  +  a2))1  ((vf+  b2))‘ 


4.132 


since  vt  a  =  v2  b  =  Vj  b  =  v2  a  =  ab  =  0 


Section  5  The  Hot-Wire  in  Compressible  Flow 


We  shall  see  that  hot-wire  fluctuation  measurements  require  detailed 
knowledge  of  the  steady-state  heat  loss  laws.  Wire  response  to  mean  flow 
is  well  defined  for  the  incompressible  case.  For  isothermal, 
incompressible  flow,  a  hot-wire  responds  only  to  velocity  changes  and  the 
output  can  be  correlated  quite  well  over  a  wide  range  of  Reynolds 
numbers.  Fig.  5.1  is  a  summary  plot  of  heat  transfer  measurements  for 
circular  cylinders  in  subsonic,  continuum  flow. 

However,  at  high  speeds,  wire  response  is  more  complex  since  wire 
recovery  factor  is  a  function  of  both  Mach  number  and  Knudsen  number. 
Fig.  5.2  was  prepared  as  a  guide  to  the  experimental  variation  of  Nusselt 
number  as  a  function  of  Reynolds  number  and  Mach  number.  The  sensor 
output  is  reasonably  well  behaved  for  supersonic  Mach  numbers  as 
indicated  by  the  lower  curve  of  fig.  5.2.  However,  the  output  is  Mach 
number  dependent  in  the  transonic  range  particularly  at  low  Reynolds 
number.  It  will  be  seen  that  the  slope  of  the  Nu  vs.  Re  relationship  is  of 
particular  concern  in  turbulence  measurements.  Fig.  5.3  shows  the 
measured  exponents  as  a  function  of  Reynolds  and  Mach  numbers  for 
several  investigations.  At  high  Mach  numbers,  the  exponent  is  seen  to 
vary  monotonically  between  the  free  molecular  and  continuum  values. 
For  an  insulated  wire,  the  slope  begins  to  deviate  from  the  continuum 
value  at  wire  Reynolds  numbers  below  200.  In  continuum  flow,  wire 
recovery  temperature  is  a  function  of  Mach  number  since  there  is  a 
changing  relationship  between  frictional  and  compression  effects.  But,  as 


,1 

A 
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Mach  number  increases,  these  effects  cancel  so  that  the  recovery 
temperature  ratio  becomes  approximately  constant  at  supersonic  Mach 
numbers  (fig.  5.4).  In  the  transitional  regime,  measurements  indicate 

that  at  Knudsen  numbers  of  about  0.1  the  recovery  temperature  begins  to 
rise  above  the  high  Reynolds  number  value.  Thus  the  recovery 
temperature  can  range  from  below  to  above  total  temperature. 
Measurements  have  been  made  over  the  complete  range  range  from 
continuum  to  free  molecular  Knudsen  numbers.  These  results  are 

summarized  in  fig.  5.5.  The  direct  effect  of  wire  Reynolds  number  on 

wire  recovery  temperature  in  supersonic  flow  can  be  determined  from 

fig.  5.6. 
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Fig.  5.1  Summary  of  heat  loss  from  circular  cylinders  in  cross-flow 
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Fig.  5.2  Empirical  correlations  of  hot-wire  heat  transfer 
at  low  Reynolds  number. 


Reynolds  number,  Re0 

Fig.  5.3  Slope  of  Nusselt  number  -  Reynolds  number  relation 


0 


A 


Mach  Number,  M 


Fig.  5.4  Recovery  temperature  variations. 


Fig.  5.5  Normalized  variation  of  recovery  temperature 
with  Knudsen  number. 


Fig.  5.6  Variations  of  recovery  temperature  in  supersonic  flow  (ref.  1). 
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The  derivation  of  the  general  fluctuation  sensitivities  of  a  hot-wire 
anemometer  involves  the  perturbation  of  the  steady-state  heat  transfer 
law  and  expressing  the  result  in  measurable  electrical  and  fluid  flow 
properties.  In  compressible  flow,  dimensional  analysis  and  test  data 
imply  that 

Nu  =  Xd/k  =  f(Re,  M,  0)  9  =  Tw/Tt  5.1 

and,  assuming  that  the  electrical  energy  input  to  the  wire  equals  the  heat 
loss  due  to  convection,  we  may  write 

W  =  I2R  =  Nut7tktL(Tw  -  r|Tt)  5.2 

Consider  first  the  electrical  properties  where 

dW/W  =  2  d  In  I  +  d  In  Rw  5.3 

But,  since  the  wire  current  does  change  with  wire  resistance,  we,  in  fact. 


dW  2  d  In  I  ^  n  .  j  «  c  / 

Tir  =  ^  p  in  Rw  +  a  In  Rw  5  4 

W  a  in  Rw 

Now  if  we  define  the  finite  circuit  parameter  e  =  (d  In  I)/(d  In  Rw)  we 
obtain 

din  W  =  (l -2e)  d  In  Rw  5.5 

or 

d  In  W  =  (1  -  2e)(d  In  Rw/d  In  Tw)d  In  Tw  =  (1  -  2e)K  d  In  Tw  5.6 
where  K  =  d  In  Rw/d  In  Tw.  Now 

d  In  e  =  d  In  I  Rw 

**  d  In  I  +  d  In  Rw 

*0 -IffJ-o-'*  5.7 

Since 

dine  -  (1  din  Tw  =  (1  -  e)K  din  Tw  5.6 

a  inTw 

then 

dine  -  (1-I-gidlnW  59 

( 1  -  2e) 

Also,  combining  equations  5.6  and  5.9,  we  determine  that 

d  In  Tw  =  d  In  e  /K(l  -  e)  5.10 

These  are  important  relationships  between  wire  heat  transfer  and 
measurable  wire  properties,  which  will  be  used  later. 

Now  we  need  to  determine  similar  logarithmic  variations  in  heat  loss.  To 
do  this  we  write 

H  =  7tLkt(Tw  -  nTl)Nu,Rel  5.11 
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Now,  using  the  chain  rule 

dH  =  dkt4r  + dTw|^  +  dNu|^- +  dn|^  +  dTt|^  5.12 
T3k*  dTw  3Nu  3n  3Tt 

If  we  divide  by  H,  the  first  term  on  the  RHS  of  equation  5.12  becomes 


J.3H  muiw-  nit;  ^ ^ 

H3ktdkt"  ktNu(Tw-  nTt)  dk*  =  kt  "  d  kt  cn<dlnTt 


Nu(Tw-  nTt) 


dk< 


the  second  term  becomes 


ktNu 


ktNu(Tw-  nTt) 


dTw 


dTw 


(Tw-nTt)  (Tw/Tt-n) 


Tw/Tt  d  InTw  =  — —  d  InTw 


(0  -n) 


¥5 

rfw. 


the  third  term  becomes 
kt  (Tw-nTt) 


ktNu(Tw-  nTt) 
the  fourth  term  becomes 


d  Nu  - 


d  Nu 
Nu 


=  d  In  Nu 


TtktNu 


k*Nu  (Tw-nTt) 


dn 


dn 


n 


(Tw/Tt-n)  (0  -  n) 


d  In  n 


and  finally  the  fifth  term  becomes 


n  k*Nu 


ktNu(Tw-  nTt) 


dTt  =  - 


n 


d  In  Tt 


d  In  Tt 


(0  -n) 

(0  -n-0) 

=  "l  (0-n)  / 

■-(1*(9^))dlnT’ 


Now  defining  the  overheat  parameter  as 

Tw  -  Tr  (0  -  n) 


'twr  = 


Tr 


n 


5-13 


and  collecting  terms  we  obtain 


I 


d  In  H  -  — d  InTw  ♦  fn»  -  1  -  Id  In  Tt 
e  -  n  1  twr  J 

-  -i-  d  In  n  +  d  In  Nu 

"CwT 


5 .14 


Now,  compressible  flow  measurements  suggest  that  Nu  =  f(Ret,  M,  0,  p,  <t>) 
and  rj  =  f (  Ret,  M,  0,  p,  <}>).  Thus,  as  a  first  step,  we  must  determine  the 
auxiliary  dependence  of  M  and  Ret  on  the  sensing  variables.  Accordingly, 


we  may  write  Ret  =  pud/pt  as 

d  In  Ret  =  d  In  p  +  d  In  u  -  d  In 
if  m  s  Ttm*  then 

d  In  Re{  =  d  In  p  +  d  In  u  -  m(  d  In  Tt 

Now  the  iso-energetic  compressible  relation  gives 
„2  u2  y-1  u2  T1  u2 


-Wf 


5.15 


5.16 
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so  that 

d  In  M  =  a_1(d  In  u  -  l/2(d  In  Tt)) 
where  a  =  ( 1  +  ^  M2  )'* 

Now,  using  the  chain  rule  once  again,  we  may  write 

.  3n  3n  3n  .  3n.A 

dn'ar4  p  ^ 

dividing  by  t],  we  may  write  the  first  term  as 

fil 1B21  =  p.n  P-  d  in  Re, 
In'  Re*  /  Re*  o  In  Re* 

and,  treating  successive  terms,  we  obtain 


5-19 


d  in  n  = 


3  Inn  ^  t,.  3  In  n  ^  ,  3lnnJ<_ 

d  In  Re*  +  ^  m  fd  In  M  +  —  d  ln  p 


3  In  Re* 

3  In  n 
3  In  <t> 


3  In  M 


3  In 


d  In  <t> 


which  may  be  rewritten,  using  equations  5.16  and  5.18,  as 
3  In  n 


d  In  n  -  ^  (d  In  u  +  d  In  p  -  m*  d  In  Tt ) 

3  In  Re*  r  T 

3  inn  /d  ln  u  dlnTt\  3  Inn 


3  in  M  '  a 
3  in  n 


2oc 


3  In  p 


d  In  p 


3ln$ 


d  In  $ 


5.20 


In  a  similar  manner  we  determine  that 

. ,  .T  3  InNu.,  _  3  InNu.,  3lnNuJf 

d  InNu  =  - — - — d  In  Re*  +  - d  InM  +  -r - d  In  p 

'  j  i  —  t  r  J  _  V 


3  In  Ret 

+  3lnNu 
+  3ln6 


3  InM 

^  ^  3  InNu  j  v 

d  In  9  +  -r -  d  In  0 

3  In* 


3  In  p 


5-2  1 


Now,  since 


d  In  9  =  d  In  Tw  -  d  In  Tt 


(e-£) 


We  may  rewrite  equation  5.21  as 
3  InNu 


d  ln  Nu  = 


- (d  In  u  +  d  In  p  -  m*  d  In  Tt ) 

3  In  Ret  r  T 

+  3hiNu  (d  ln  Tw  .  d  ln  Tt )  +  d  In  p 


3  In  9 

3  InNu 
3  ln  $ 


3  ln 


d  ln  $ 


5-22 


5-23 


Substituting  for  d  ln  r|  and  d  ln  Nu  in  equation  5.14  and  collecting  like 
terms,  we  determine  that: 
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,  .  ..  f  Ret  3Nut  0  3  NU* 

d  In  H  -  n*  -  mt  — -T  ■  -  - - 5-— 1 

1  *  1  Nut  3  Ret  Nut  3  6 

1  r«  1  M  3n  _  Ret  3n  'i  1  M  3Nuti  J T 

■  “J1 '  t  m ' m<  ~  7r^->  '  2^Fu,im  j d  ln  T' 

fRet3Nut  1  M  3Nut  1  r  1  M  3n  Ret  3n  ^iJU 
+  l  Nut  3  Ret  +  "oT  Nut3M  "  t^A  a  n  +  n  3ReA-ld  nu 

♦  f  Bfj^Q-ldlnp  J  ®- 4"H»._®_1dtalV 

l-Nut3Ret  twr  n  3RetJ  r  I  Nut  3  0  0-rw 


[♦  3Nut  1  ♦  3n  1 
Nut  ~3$  rT  1H  J  d  ln  * 


Since  we  have  previously  shown  that 
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d  ln  W  =  d  ln  H  -  (V-^l  d  In  e  and  d  ln  Tw  =  S~--; 

\  l  -  e/  K(l-c) 


we  can  rewrite  equation  5.24  as 


(Ll23.\  d  in  e  -  r^LiS-Nu  +  _01  dine 
vi-e2aine  Ldine  e-nJKd-c) 

r  Ret  3  Nut  0  3Nut 

Lnt  m^Nut3Ret  Nut  3  9 

1  f,  1  M  3ti  ^  Ret  3n  ^  1  M  3Nu^i  ^ ^ 

~  t^A  "  2c c  n"  3M  ’m,7  3  Reti  '  2a  Nut3M  J  d  n  * 

fRet3Nut  1  M  3Nut  1  fl  M  3n  Ret  3r>  ^  J(_ 

+  I  Nut  3Ret +  !TNut3M  n  3M  +  n  lReAJdlnU 

rRe,  |Nut_J_  Re,  3n_ldln 
l  Nut  3  Ret  twr  n  dRe,J 


Now 


r  ♦  3Nut  1  4>  3  n  1 

twr  n  3<t>  J  ln  ^ 


d  ln  W  -  d  ln  I2Rw  -  d  ln  I  +  2(d  In  I) 

3  ln  I 
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and,  for  fixed  flow  conditions,  equation  5.14  gives 


d  In  H  =  [ 


d  In  Nu  _ 
d  In  e  T  0 


e 


hi 


d  In  Tw 


5.27 


So  that,  introducing  Aw'  =  1/2(9  In  Rw/9  In  I)  and  equating  W  = 
obtain 


(2A'W  +  2)  din  I 


fd  In  Nu 
L  din© 


e 

e  -n 


] 


2A; 


K 


w  d  In  I 
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so  that 


d  In  Nu  _ 
d  In  0  T  0 


e 


hi 


5.29 


and  so  the  LHS  of  equation  5.25  reduces  to 


ln*-K[w^] 


d  In  e 
K(1  -6) 


If  we  define  E'  =  (1  -  e)/(l  -  2eAw’),  then  we  determine  that 


d  In  e  «  E'A; 

_1_ 

twr 


v  [[  nt  ~ 


m« 


Re,  3Nu,  0  3Nu, 

*  Nu,  ^Re,"  Nu,  a  e 

f  1  M  3n  Re,  3n  ^  1  M  3Nu,i 

l  2an  aM  t  n  3Re,J  ~  2aNu,3M  Jd  * 


r  Re,  aNu,  1  M  3Nu,  1  /  I  M  3d  Re,  3n 
l  Nu,  3  Re,  +  a  Nu,3M  ~  t^llTT  9M  +  ~  3Re,;Jd  nU 


rRe,3Nut  i_  Rfj  3n_idln 

Nu,  3  Re,  twr  n  3Re,J 

r  ♦  3Nu,  1  4>  3n  i  .  1 
^Nu,  dt  twr  n  3$  ■*  n<M 
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where  the  last  term  in  equation  5.30  may  be  rewritten  as 

fj_  3.1nn  .  3lnNu,ld 

*■  twr  a  4>  a  4>  J 


H,  we 


Now,  following  Kovasznay  (ref.  9),  the  basic  equation  for  a  hot-wire 


531 


inclined  to  the  flow  may  be  written  as 


e' 

e 


P  u  Tt 
Ae«  —  +  Ae„—  +  Ae_  ±  Ae„ 

P  p  “u  TTTt  4 


v_' 

u 


where  the  sign  conventions  must  be  determined  separately  for  constant 
current  and  constant  temperature  operation  as  will  be  shown  in  section  6. 
Replacing  the  logarithmic  variables  by  fractional  perturbations  and 
approximating  d<|>  by  v'/u,  the  hot-wire  sensitivities  in  equation  5.31  may 
be  finally  written  as 


Ae 


3  In  Nu\  ( 
3  In  Ret/~  \ 


3 

3 


In  Nu 
In  0 


) 
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J_  fi  m  l  /3lnNuV>1 

twrV1  2oc  n  3M  *  n  3  Ret-J  2cA3lnMA>J 
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and 


Aev  = 


E' 

100 


1  f  3  In  n  ^ .  /3ln  Nu\ ' 
"Cwr  ^  3  $  -'''34*  ' 
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We  shall  see  that  the  proportionality  constants  are  related  to  the 
particular  electrical  system  and  are  different  for  constant  current  and 
constant  temperature  applications.  However,  these  hot-wire  sensitivities 
apply  to  all  flows  whether  subsonic,  supersonic,  incompressible  or 
compressible  whether  continuum  or  free  molecular.  For  supersonic  Mach 
numbers  ie.  M  sin  <}>  >  1.2  we  have  seen  that  all  derivatives  with  respect  to 
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Mach  number  are  negligible  so  that  Aep  =  Aeu.  Now,  since 

5(pu)  3p  3u 

pu  p  *  u 

equation  5.31  may  be  rewritten  as 

TV  (pu)‘  v‘ 

e  =  Ae.  ^  +  Ae_..-^rr  ±  AeM— 


'tt  Tt  p  u  *  *  u 


5.37 


Now  let  us  consider  how  to  evaluate  various  fluctuating  terms  for  the 
case  of  a  normal  hot-wire,  where  we  can  rewrite  equation  5.37  as 


TV  (pu) 

e  =  Ae_  -2-  +  A©  “zzr 
Tt  Tt  pu 
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We  start  by  taking  into  consideration  the  so  called  Kovasznay  diagram 
(ref.  9)  as  follows.  Squaring  both  sides  of  equation  5.38  we  get 


•*75  ,  TV2  (pU)'TV  2  (£U) 

6  "  AeT.  T,2  *  2  (AV  (p51)  *  AV  (pn)2 

Dividing  through  by  (AeTt)2  the  above  expression  becomes 
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e 

(A© 


2  (pu)VAepu\2  ^(pu)'TV TV2 

Tl>2’  (pu)2VAei/  +  2<P3T«)  '■^®Tt '  *  Tt2 
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If  we  now  define 

-H 

<AV 


Ae. 


2  =  s' 


and 


pu 


Ae, 


=  r 


Tt 


then  our  expression  can  be  rewritten  as 
^2 


(pu)~ 

(f?u) 


2  r  +  2  (^uTt) 


(pu)'TV  TV2 

r  +  1fT2 
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Kovasznay  then  suggests  we  plot  s2  versus  r.  This  is  a  second  order  curve 
in  r  and  the  coefficients 
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(f>u)'2  (pu)'Tt  TV2 

^=")2 '  Wnj  '  T2 

are  obtained  by  a  least  square  fit  of  the  function  s2  =  f(r).  These 
coefficients  then  become 


where 


N 

N 

S  S,2 

i  =  1 

N  , 
S  r,2 
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(pu)*TV 
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1 

D 
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1=1 
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n  , 
S  r,3 
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i  =  f 

n  ,  , 

2  r(2s,2 

1  =  1 

N  - 

S  r,4 

i  =  1 

N 

Ss,2 

i=i 

N 

S  r, 

1=1 

N  , 
S  r,2 

i  =  1 

TV2 

Tt2 

1 

D 

N 

S  r,s,2 

i  =  1 

n  , 

S  r,2 

i=l 

N  , 
S  r,3 

i=i 

N  ,  , 
S  r,2s,2 

i=i 

S  r,3 

i  =  i 

S  r,4 

i  =  1 

N 

N 

S  r, 

1=1 

N  , 
S  r,2 

i  =  1 

D 

= 

N 

S  r, 

1  =  1 

N  , 

S  r,2 

i  =  1 

N  , 
S  r, 3 

i  =  1 

N  * 

S  r,2 

i  =  1 

N  * 

S  r,3 

i  =  1 

N  * 
S  r,4 

i  =  1 

and  N  is  the  total  number  of  points  used. 

In  order  to  obtain  other  terms  which  appear  directly  in  the  turbulent 


momentum  and  energy  equations,  we  assume  that  we  have  an  isentropic 
flow  field.  This  permits  us  to  write  the  energy  equation  in  its  differential 
form  as 

1  dTt  ,  .v,,2du  dp  dp  c 

oc  Tt  U  P  P 

where  oc  ■  1/(1  +  Y  ~  *M2) 

Next  let  us  consider  the  equation  for  the  mass  flow  per  unit  area  and  time 
in  its  differential  form 


dm  3u  |  dp 
m  “  u  +  P 

Substituting  for  3p/p  in  equation  5.42  gives 

1  _  ,  ,,,,2  3u  .  3p  3m  3u 


m  +  U 


oc  Tt  UP 

or,  collecting  terms,  we  obtain 


3u _ 1 _ dTt _ 1 _  rdp 

u  a  [1  +  (v  -  1)M2]  T*  [1  +  (y  -  1)M2]  *■  P 


5.43 


5-44 


5-« 


Which,  defining  (3  =  a(y  -  1)M2,  can  be  written  as 


(  1 

-t 

'  a 

_u§p .  < 

—  = 

U 

'OC+  3 

)  Tt  ' 

'  a  + 

P 

so  that,  since 

p’/p  ^ 

1.0 

u’ 

<  1  ' 

+  (  «  'I 

G 

'oc  +  3 ' 

hi 

'  oc  +  3' 

(pu) 


5.46 
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Squaring  both  sides  of  this  equation  leads  to  the  expression  for  the 
streamwise  turbulence  intensity 


(  1  1 

f  2a  1 

.(pun,- 

r  «  f 

\oc  * 

1  TP  1 

^  a  +  B^ 

'(puTt)  1 

5-46 


In  order  to  obtain  an  expression  for  the  density  fluctuation,  we  consider 

the  mass  flow  equation  in  the  form 

3p  _  3m  _  3u  c  4Q 

p  ~  m  u  y  y 
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Substituting  in  equation  5.46  and  again  neglecting  pressure  fluctuations 
we  obtain 

3p  r  oc  1 3  m  /  I  \  3  T* 

P  L1  ~  a  +  B  J  m  ~  +  B '  Tt  5  50 

which  may  be  written  as 


f  8  ' 

v3(pu) 

-f  1 

\3Tt 

\a  +  p  > 

'  pu 

Va+  B 

'  Tt 

(  *  ) 

(pu)' 

r  1  ) 

Tt 

'oc  +  B/ 

pG 

B/ 

Tt 

while  the  squared  term  may  be  written  as 

P^2  _  /  B  \2  (pu)'2  /  2B  \(puTiV  +  /  1  \2  TV2 
p2  =  U  +  B'  (pH)2  "  'a  +  B'(puTt)  +  +  B'  Tt2 


The  mass  fluctuation  can  now  be  derived  by  setting 

Pu’  m£x  V  m  \(  P  \(Pu)‘  (  1  ^  Tt  1 

pu  P  u  a  +  B '  pu  '  a  +  B  '  W  ■* 


p  u  L\  a  +  B  ^  pu  '  a  +  B  '  Tt  ■» 
LVa  +  B^pu  \  a  +  B  '  Tt  ■* 


time  averaging  we  obtain 


ocB  (pu)'z  +  (B  -  oc)  (pu)'Tt'  /  1  \2  Tt'2  _ 

:  +  B)2  (pu)2  +  (oc.  *  B)2  (puTt)  'a  ♦  B '  Tt2 


pu  (a 


(puTt) 
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We  can  also  determine  expressions  for  the  static  temperature  fluctuation 
and  the  term  u'T'  which  is  obtained  by  considering  the  perfect  gas  law 

p  =  R(pT)  5.56 

or 


3t  3p  Bp 

T  =  p  ”  p 


5-57 


Neglecting  9p/p  because  it  is  small  compared  to  the  other  terms,  we  see 
that 


Bt  Bp 
T  =  '  f> 


T‘  p ■ 
or  =-=--£• 

T  p 
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then  from  equation  5.53  we  see  that 

T‘  .  \j  B  \  <pu)'  _  /  1  \Tt1 
T  Lxa+g/  pu  '  a  +  3  '  Tt  J 
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so  that 


UTS  U'  x  T'm  W  oc  ypu)',/  1  \Ttl 
uT  n  T  l-xa  +  0/pu  \a  +  3'Tt^ 


\(  e 

\<pu>’  1 

r  1  1 

|ll 

LVa+  3 

'  pu 

<cc  +  BJ 

1  Tt  J 
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which  may  be  rewritten  as 


u'T'  =  r  («  -  3)  (pu)'Tt'  0 eg  (pu)'z  +  1  Tt2 1  r6l 

uT  L(a  +  B)2(puTt)  “  (cc +  3)2  (pu)2  +  (cx  +  3)2  Tt2  J 


(pu)‘ 


The  relationship  between  p'u'  and  u'T’  can  be  determined  using  equation 
5.58  since 

5-62 


-£'*-u-’=  -lx*’ 
p  u  T  u 


so  that 


76 
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Now,  since 


„  u'  T 

B  +  a  = 
U  T 


(a  +  B)  I  +  8 

T  (pu) 
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so  that 


=  ii  _  X  5  70 

(pu)  G  T 

It  should  be  noted  that  this  last  expression  can  also  be  derived  from  the 
fact  that  m  =  pu  and  p  =  RpT  provided  the  assumption  is  made  that  the 
p'/p  is  negligible. 

It  is  important  to  recall  that  all  the  above  relationships  can  be  derived 
from  single  normal  hot-wire  measurements.  However,  other  turbulent 
terms  which  appear  in  the  momentum  and  energy  equations  require 
inclined  hot-wire  measurements. 


Let  us  now  consider  the  case  of  a  yawed  hot-wire.  If  we  return  to  the 
hot-wire  equation  5.37  and  square  it  we  see  that 


e2  = 


AeTT  Tt  +  Ae^  (pu) 


<£">' .  Ae..  1' 


V  U 
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expanding  and  collecting  like  terms  we  obtain 


(AO2  *  <AeT/  p  *  <A*V>2  |2  *  2 


(pu) 

* 

% 


(pu)'Tt' 

(pu)Tt 


+  2  (AeT )(Aev)^-  +  2(Ae^u)(Aev) 


(pu),v' 

(pu)u 


5.72 


Now  we  run  the  test  twice  at  each  point,  once  with  the  hot-wire  at  <J> 
and  a  second  time  with  the  wire  at  <f>  =  180°,  then  we  can  write  that 

_  n 


-  0° 


e i=Q'~  %=160e 


since 


5?3 
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•SI 
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•3 


1  ’i.l 


rt 

a? 


a 


»!*? 


$ 


(Avtee-)S' ]  •  K^S' ] 


0=180° 


(A6t )(Aev)^-|  ■=  -  |(Ae  >(Ae„)?J| 

T'  uTtJ».=o*  1  '  uT*J*=i»o- 


[(AVKAe.)^]  - 

I  (pU)U  Jtf_QO 


(8%)(*ey,<£HIr 

P  y  (pu)u  J^.,80. 


Dividing  through  by 


±  4  (AelT)(Aev) 


we  obtain 


a,‘2  _  p'2 


0=0°  g0=iaoa 


±  4  (AeT  )(Aev) 


f  A%yp^)'v'  +  W 

=  \Ae,  Apm  +  uTt 


However,  if  we  define 


i  2  -  ?2 

e0=O°  %=  180° 


±  4(Ae-)(A©„) 


then  equation  5.74  can  be  rewritten  as  follows 


s*  =  <£H>2'  r  *  ?£ 
(pu)u  uTt 
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If  we  now  plot  s*  versus  r  we  see  that  the  function  is  linear  with  a  slope 


•  »,<.)  ‘ 


) 

and  an  intercept  of 


Tjv' 

uTt 


These  values  can  also  be  obtained  if  we  perform  a  least  square  fit  to  the 
data.  It  can  be  shown  that,  if  we  have  N  number  of  data  points,  this 
method  gives 


(pu)'v' 

(pu)u 


N  , 

NZn  s?2 

4  =  1 


N 


N 


X  r<  X  S  r 


1  = 1 


181 


N 


N  X 
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and 
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At  this  point  it  should  be  noted  that  from  the  normal  hot-wire 
measurements  we  have  values  for 

Pf,  .  and  £ 

(pu)  (pu)Tt  Tt 2 


and  from  the  yawed  wire  measurements  we  now  know 

tegZ  and 

(pu)u  nTt 


2 17.2 


Since  e'^.Q  will  also  have  been  recorded  we  can  obtain  the  value  of  v  /u 


from  the  following  expression 


We  can  now  determine  the  Reynolds  shear  stress  term  u'v’  by 
multiplying  equation  5.47  by  v'/u  to  obtain: 


u'v 

-  (  1  'i'K7  . 

(  «  \(pu>-v 
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G2 

^a+p/GTt 

\0C  +  B'  (pG)G 

The  term  v'T' 

which  appears 

in  the  turbulent 

Prandtl  number 

formulation,  equation  3.22,  is  obtained  by  multiplying  equation  5.59  by 
v’/u  to  obtain: 


v  T  =  -  1 

/  B  \  (pu)'v' 

c  1  'I  T<v'  1 

GT  1 

■'a  +  B'  (pu)u 

''a  +  B 1  uTt  J 

Two  other  useful  relationships  may  also  be  derived. 
5.67  by  v'/u  gives 


Tiv' 

GTt 


=  B 


u  v’ 
G2 


TV 
a  — =■ 
GT 


Multiplying  equation 

5-63 


Similarly,  from  equation  5.70  we  obtain 

(pu)'v'  _  uV  _  vT 
(pu)u  G2  GT 
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The  triple  correlation  p’u'v'  which  appears  in  the  compressible  shear 
stress  formulation  (see  section  1)  can  be  determined  from  equation  5.47, 
5.52  and  the  rms  measurement  of  equation  5.80.  So,  at  this  point  we 
know,  in  principle,  how  to  obtain  the  turbulent  flow  properties  and 
correlations  which  appear  in  the  momentum  and  energy  equations  from 
variables  determined  directly  by  the  hot-wire. 
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Although  many  years  of  effort  have  been  expended  in  hot-wire 
anemometry  research,  it  is  still  clearly  an  inexact  science.  We  have  seen 
that  numerous  assumptions  must  be  made  to  estimate  the  fluctuating 
flow  variables  from  the  measured  hot-wire  quantities.  Now  let  us  turn 
our  attention  to  some  of  the  practical  aspects  of  hot-wire  anemometry. 


Section  6  Practical  Procedures  for  Turbulence 
Measurements 


We  have  seen  in  section  4  that,  with  appropriate  electrical  compensation, 
the  time  constant  can  be  restored  to  frequencies  sufficient  to  assume 
essentially  instantaneous  response.  Imposing  this  constraint,  we  may 
rewrite  the  general  expressions  for  wire  sensitivity  (ref.  1)  for  constant 
current  as 


Ae  =  -  &ep  | 

o 

o 

I  -  Aeu  | 
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|  +  Ae^l 
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with 


t&p 


Ae< 


eE' 

100 


Aw 


9  In  NUo 

9  In  Re0 


Ae 


T. 


‘  ^5 [A"{ 

Aw  /  1 
twr 

-  E  [K^{K-I-n,-m‘ 


Aw  9  in  n  1 
twr  9  In  Re0  i 

■} 

.}] 


9  In  Nup  +  _1_  9  in  Nu 
9  In  R©o  a  9  In  M 


9  In  n  9  In  n 


9  In  M  9  In  Re, 


9lnNu0  1  9lnNu 


+  — 


9  In  Re0  2 a  9  In  M 


■) 


Aw  (  J_  9ln_n_+  9  inn  \1 
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6.4 
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©  E  Aw 
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where  E'  =  (1  -  e)/(l  +  2Aw'e) 
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For  constant  temperature  systems,  the  sensitivities  are  usually  defined 
with  respect  to  a  fluctuating  current,  so  we  may  write 

Al  =  Al^(lOO^)  +  AIb(iOO^t) 

-  AlT#(lOO^-°)  +  Al*(l00-sr)  6.7 

We  can  obtain  the  sensitivities  for  A I  from  those  for  Ae  by  replacing  eE' 
by  IE”,  where  E"  is  the  feed-back  reduction  factor  defined  in  (ref.  10)  as 

E" - l—r  6.6 

2  Aw  +  j 

Thus,  for  the  correct  interpretation  of  the  electrical  signal  from  a 
hot-wire,  we  must  consider  the  response  of  the  electrical  system,  the 
variation  of  wire  resistance  with  operating  temperature  and  the  fluid 
dynamic  calibration. 

Let  us  first  consider  the  electrical  system.  We  can  see  that  an  important 
term  in  any  hot-wire  anemometer  measurement  is  the  feedback 
parameter  (e)  which  relates  wire  power  loss  to  measured  wire  voltage. 
Consider  then  the  constant-current  circuit  shown  on  the  next  page  where 
e  =  -(d  log  Iw)/(d  log  Rw)  and  its  equivalent  circuit  derived  using 
Thevenin's  theorem. 


Now  the  open-circuit  voltage  across  AB  is 


V (100  +  Rb) 
(Ri  +  Rb  +  100) 


6.9 


and  the  resistance  in  series  with  AB  when  V  is  shorted  is 

ro.  ( — L  +  JLy1 

U00  +  Rb  Ri) 

m  Ri(Rb  +  100) 

(Ri+Rb+100) 

so  that  the  wire  current  is 

_  Vo  _ V  (100  +  Rb)(Rh-Rb  +  100) _ 

iw  Ro  +  Rw  +  10  =  (Ri+Rb+100)(Ri(Rb+100)+(10+Rw)(Ri+Rb+100)) 

_ V  (100  +  Rb) _ 

Ri(Rb  +  Rw  +  1  10)  + (10+ Rw)(Rb+ 100)  011 


Since 

d  log  Iw  -  Rw  (Ri+Rb+100) 

d  log  Rw  =  Ri(Rb  +Rw  +  1  10)  +  (10+Rw)(Rb+100) 


then 


_ Rw  (Ri+Rb+100) _  ,  i 

€  "  Ri  (Rb  +  Rw  + 1 10)  +  (10+  Rw)(Rb+  100)  '  3 

Now  consider  the  constant  temperature  system  shown  in  the  figure  below 
where  the  amplifier  is  characterized  by  its  transconductance,  G,  and  its 
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short  circuit  current  Isc. 


Now  I'  =  GAe  +  Isc  and  the  relationship  between  the  current  I'  supplied  to 
the  bridge  and  wire  current  Iw,  is 


Iw 


roc’  where 


R2+  R3 

Ri+  R2  +  R3  +  Rw 


6.14 


The  unbalanced  voltage  is  proportional  to  the  bridge  unbalance  Rb  -  Rw  so 
that 


Ae  =  <xT(Rb  -  Rw) 


6.15 


Perturbing  these  equations  leads  to 

3  log  I'  3  log  Iw  R wGoc' 

e  =  -  — -  =  -  — -  =  - - - -  6. 16 

3  log  Rw  3  log  Rw  1  -  Goc'(Rb  -  Rw) 

Thus  we  can  relate  wire  voltage  to  the  wire  input  power  for  both  constant 
current  and  constant  temperature  anemometer  systems.  Next  we  need  to 
determine  the  relationship  between  the  measured  output  voltage  and 
actual  wire  voltage  as  shown  on  the  next  page 
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Now  let  us  look  at  the  parameters  which  govern  the  relationship  between 
wire  resistance  and  operating  temperature.  First  consider  the  resistance- 
temperature  relation  and  some  wire  properties  given  below. 
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Wire  Properties 
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The  resistance-temperature  relationship  may  be  expressed  as 

Rw  =  Rf{l  +  af(Tw  -  Tf)  +  yf[af(Tw  -  Tf)]2}  6.21 

where  f  is  some  reference  condition.  The  parameter  K  which  relates  wire 

power  input  and  temperature  can  be  calculated  from 

„  1  /  ^  x  d  log  Rw 

K  =  —  (ar  Tr  +  aw)  =  .  .  _  T 

1  +  aw  d  log  Tw 

Rf 


=  —  Tw  [at  +  2Yfccf2(Tw-TO] 


~  oct  Tw  [  1  +  2  y  toct  (Tw  -  Tf )  ] 


6.22 


Plots  of  K  and  wire  temperature  versus  wire  resistance  ratios  are  shown 
in  figures  6.1  and  6.2. 

Wire  overheat  parameters  may  be  determined  as  follows 

aw'  =  (Rw-  Rr)/Rr  6.23 

where  Rr  is  the  recovery  resistance  of  the  wire.  The  recovery  factor 
Tr/Tw  varies  between  0.96  and  1.1,  but  is  essentially  constant  (0.96)  for 
Mach  numbers  greater  than  2.0  and  Reynolds  numbers  greater  than  20. 
Typically,  Aw'  >  aw'  by  up  to  20%  at  high  overheat  ratios  and  there  are 
several  ways  to  compute  its  value.  A  method  which  eliminates  the 
differentiation  of  experimental  data  with  a  steep  slope  is  detailed  below. 
Experimentally  it  is  easier  to  determine  Aw'  beginning  with  a  plot  of 

T(L  Rw  _  «.  Rw  -  Rr 

.  ft*,  vs  aw  =  — = - 

aw  Rr 

Now,  defining  C  as  the  slope  of  (Iw2Rw)/aw'  vs  aw'  at  high  overheat, 
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V? 


ie.  C  = 


/Tw  Rw'j  _  Iv  Rw 

\  aw  /  aw  =  0 _ aw 


6.24 


then 


Aw 


aw 


1  -  ( 1  -  aw) 


Caw 


6.25 


1  -  Caw 

We  are  now  in  a  position  to  determine  the  finite  circuit  parameter  E'. 
remaining  terms  in  the  sensitivity  equations  are  defined  as  follows 
e~  =  mean  voltage  across  the  wire  H  =  total  power  to  the  wire 

6.26 


-  aVvr) 

Tr  arTr 


cc  - 


1 


1  +  M2 


6.27 


m<  «  |log  M  =  .65  and  nt  =  =  .75 

d  log  Tt  3  log  Tt 


Fig.  6.1  Variation  of  parameter  K  with  wire  temperature 
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Fig.  6.2  Variation  of  wire  temperature  with  wire  resistance 


For  constant  current  operation  we  also  need  to  determine  the  wire  time 


constant  from  ref.  1  as 


M  * 


Cpwpwd^  Rw 
4  NUomkoair  Rr 


6.26 


A  plot  of  NuQ  versus  Reo  is  given  below.  But  we  should  use  measured  values 
of  Nusselt  number  (Nuom).  Typical  comparisons  between  measured  and 
theoretical  values  are  shown  in  figure  6.3.  Figure  6.4  shows  the  effects  of 
the  wire  dimensions  and  thermal  properties.  If  we  know  Cpw,  pw,  d2,  Nuom, 

ko,  and  Rw/Rr,  we  may  calculate  the  true  time-constant  and  correct  the 
measured  wire  voltage  as  follows 


Ae2  ( Ae2  -  Ae2^ 
Ae^  VAe2 -  Ae2 


6.29 


where  Aeo  is  uncompensated  wire  voltage  which  is  small  for  usual  values  of 


iT.VlvVi.'i.vVAVj, 


\i.  Thus, 


Ae 


Ae 


m 


„  A 
Mget 


where  psel  is  the  compensating  amplifier  time  constant.  So  we  are  now  in 
a  position  to  account  for  the  effects  of  the  electrical  system  and  wire 


resistance  variations  on  the  measured  electrical  signal. 
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Now  let  us  direct  our  attention  to  the  remaining  terms  in  the  sensitivity 
equations.  Following  the  results  of  section  5,  we  may  make  the  following 
high  Mach  number  approximations  for  M  sin  <{>  >  1.2  and  (3  In  Nuq/3  In  M) 


3  In  Nile  _ 
3  In  R©0 


/  Re0*m 
Re0  sin  ^ 


6.30 


where 


with 


Re0*rn  -  [1.1  -  .62^(/Re;  -  1.1)  5  ] 2 


u,  _  /ST 

L  V  ko,ir 


t|  =  0.95  +  0.0963/Ren 


3n 

3Re0 


.0963 


3  in  n 
3  In  Re0 


1+10  Re0 


6.33 


6.34 


Normally  we  assume  that  for  ReQ  >  40 

3  In  n  =  Q 
3  In  Re0 


6.35 


Finally  for  M  sin  <j>  >  1.2 


3  In  n 
3ln  M 


so  that  for  high  ie.  sin  <J>>1.2 


Ae  =  Ae  =  Ae 

p  u  pu 


and  we  may  rewrite  the  sensitivity  equations  as 


6.36 


eE’  r  d  m  Nu0  Aw  d  in  n  1 

Ae'u  “  TOO  LAV  3  in  Re0  "twr  3  In  Re0J 


3  In  Nu0  Aw  3  In  n 
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a  In  Nu, 
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Ae 


eE'Aw  r  3  In  H 
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100 


r  d  in  H  1 


6.36 


Thus,  for  M  sin  $  >  1.2  and  Re0  >  40,  it  is  possible  to  calibrate  a  wire 

directly  in  a  flow  with  varying  mean  mass  flow  distribution  such  as  a 

boundary  layer  or  by  varying  the  freestream  total  pressure.  With 
pud  .  I2  Rw  I2RwRr(Xr 


Reo  1 1 


and  Nun  = 


Jtlkx(Tw-Tr)  Tllk^Rw-Rr) 


6.39 


the  term 
9  lnNu0  _  Re0  9Nu0 
d  In  Re„ 


Nu0  3  Re0 
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[  2  Rw  Rr  /  (  R w  -  Rr 
So  that  if  we  measure 


-x 


ai2  RwRr/kT(Rw  -  Rr) 
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6.40 


I2  RwRr 


kT(RW  -  Rr) 


at  several  different  values  of  pu/p  and  differentiate  the  data  we  can 
determine 


a  In  Nu, 


3  In  Re, 


directly.  It  helps  to  plot 


I2  RwRr 


kT(Rw  -  Rr) 
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against  the  square  root  of  pu/p  since,  for  Re  >  40,  this  plot  is  usually 
linear  and  may  be  smoothed  before  differentiation. 


Other  forms  of  the  sensitivity  equations  can  be  introduced  which  also 
allow  for  direct,  in-situ  calibration.  If  we  put  a  normal  wire  into  a 
uniform  free-stream,  equation  6.1  reduces  to 


e' 


6.41 


Similarly,  holding  other  terms  constant,  we  can  show  that 


Ae 


pu 


and 


Ae.. 


write 
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In  pu 
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6.42 


6.43 


6.44 


6.45 


The  sensitivity  Aepu  can  then  be  evaluated  from  the  wind  tunnel 
calibration  of  a  wire  on  the  tunnel  center-line  by  varying  po  for  constant 
M  and  To.  We  also  must  hold  the  overheat  constant  at  various  values.  If 
we  record  values  of  E  versus  po  at  various  constant  overheats  and  plot 

In  E  versus  In  po  as  shown  on  the  next  page, 


It  is  recommended  that  Ae0  is  obtained  by  changing  wire  angle  at 

constant  overheat  in  a  fixed  free-stream  flow.  We  note  that  this 

calibration  should  be  conducted  at  a  variety  of  tunnel  total  pressures 

with  Mach  number  and  total  temperature  held  constant.  The  sensitivity 
varies  non-linearly  with  Reynolds  number  but  is  independent  of  Mach 
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number.  A  typical  calibration  is  shown  below. 


The  sensitivity  AeTl  presents  a  slight  problem  since  it  cannot  normally  be 
determined  in-situ  using  the  simplified  equation  6.41  as  it  requires  a 
tunnel  center-line  calibration  in  which  M  and  po  are  held  constant  as  total 
temperature  is  changed.  This  is  generally  not  feasible,  since  changing  To 
will  usually  change  po  as  well  as  pTu.  However,  for  M  sin  0  >  1.2  and  Re0  > 
40  we  can  rewrite  equation  6.37  as 

AeT_  -  ^[k*AV<K-1.75>  -  6.46 

so  that  the  total  temperature  sensitivity  can  be  calculated  using  the 
in-situ  determination  of  mass  flux  sensitivity  and  equation  6.22. 

Once  the  hot-wire  has  been  calibrated,  the  usual  test  procedure  is  as 
follows.  First  we  record  Iw,  Vw,  V'out  (compensated),  V’baiancc 

(uncompensated),  Rb,  Ri  (setting  for  overheat),  Rw  cold  and  u  .  for  each 

SCI 
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overheat  setting.  Next  we  plot  Rw  vs  Iw  and  calculate  A’w  vs  Rw.  Then,  for 
each  Rw,  we  calculate  the  following  terms 


1.  e 

.  3  In  n 
°'  3  In  Re 

15.  A 

2.  Aw 

9.  e  -  wire  voltage  (Vw)  16.  Vo«t 
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17  Vout 
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4.  K 

1  1.  Aen\ 

16.  e' 

5-  3lw 

12.  AeT# 

19.  elsorr 

6.  "twr 

13-  Iw 

2  0.  r  the  sensitivity  ratio 

j  3lnNu 

3  In  Re 

14.  Cw 

Aem 

AeT, 

If  we  now  take  the  root-mean -square  hot-wire  signal  and  divide  by  AeTo 
we  obtain 

+  r2(p  u)’2  -  2r^%  6.47 

and  we  may  determine  the  three  unknowns 

T?,(pu)'2  and  (p u)'T; 

In  principle,  only  three  overheat  settings  are  required.  But,  in  practice,  at 
least  three  times  that  number  should  be  made  in  order  to  provide 
consistency  checks  and  reduce  the  scatter.  The  results  can  also  be 
determined  from  the  so-called  Kovasznay  diagram  (ref.  9)  or  by 
regression.  The  Kovasznay  diagram  can  be  obtained  by  plotting  o  = 
e’/AeTo  against  the  sensitivity  ratio  (r).  Typical  curves  showing  the 

dependence  on  the  relative  magnitude  of  the  three  unknown  quantities 
are  presented  in  figure  6.5.  If  only  mass  flow  fluctuations  are  present 
the  diagram  will  be  a  straight  line  from  the  origin  with  the  slope 


proportional  to  the  mass  flow  fluctuation  level.  If  only  total  temperature 
fluctuations  are  present,  the  diagram  is  a  horizonal  line  at  the  total 
temperature  fluctuation  level.  If  both  fluctuations  are  present,  the  plot 
varies  according  to  the  correlation  coefficient  between  the  two 
fluctuations.  For  perfectly  correlated  or  anti-correlated  fluctuations,  the 
curve  becomes  a  straight  line,  and  when  there  is  zero  correlation  it  is  a 
hyperbola.  In  these  cases,  the  two  fluctuation  levels  can  be  determined 
as  indicated  in  figure  6.5. 


For  a  yawed  wire  we  have 


e'  =  Ae  ~  +  AeT  -5-  ±  Ae.  — 

pu  pu  Tl  Tt  ♦  u 


6.48 


so  that,  when  we  difference  the  mean  square  of  two  readings  taken  180 
deg.  apart,  we  obtain 


e’2=  e'02o  -  c,280.  =  4[AepuAei(pu)’v'  +  A  e0A  eTtv'Tt’]  6.49 


which  we  may  write  as 


s*  = 


4  A  e  ^  A  e  Tt 


=  4  [v’Ti  +r*(pu)’v’] 


6.50 


where 


_  Aep-Ae» 


r*  = 


A  e  0  A  e  Tt 


6.51 


If  we  plot  s*  against  r*  we  obtain 


(pu)’v’  and  v’T,’ 


from  the  intercept  and  slope  respectively. 


Now,  having  obtained  the  rms  values 
(pu)'  Tt 


,  -=-  and  the  correlation  (pu)'Tt’ 
pu  Tt 


from  a  normal  wire,  and 

(pu)'v'  and  VTt 

from  a  yawed  wire,  we  are  in  a  position  to  determine  the  fluctuating  flow 
variables  outlined  in  section  5. 


Unfortunately,  the  mass  flow  and  total  temperature  fluctuations  are  not 
the  sole  characteristic  parameters  of  an  unsteady  supersonic  flow  field. 
We  must  also  consider  the  vorticity  t,  entropy  a,  and  sound  %  modes. 


If  the  fluctuation  levels  are  low,  as  in  the  free  stream,  these  modes 
satisfy  separate  linear  differential  equations  and  the  hot-wire  equation 
becomes  (ref.  1) 


where 


Ae  =  a  Ae„  +  t  Ae,  +  n  Ae„ 

o  x  n 


6.52 


Ae„  =  Ae„  +  a  Ae„ 

C  p  T 

AeT  =  (3  Aej  -  Aeu  =  p  AeT  -  Aep  (M  >  1 .2) 
Ae„  =  cc(y  -  1)(1  -  nxM)  AeT  -  ^  Ae„-  Aep 

s  oc(y  -  1)(1  -nxM)  A©t  -  (j*  +  l)Ae^  (M  >  1.2) 


6.53 

6.54 


6.55 


a=i/<i  +  r-T 


M2)  =  (1  +  0.2M2) 


where 


B  =  a(v  -  1)M2  -  0.4aM 


for  air  and  nx  is  the  direction  cosine  defined  as  nx  =  -1/M. 


The  fluctuation  variables  are  given  in  ref.  1  as 

=  <r  +  (y  -  l)ir 
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=  vrr 


1  x 
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6.56 

6.57 

6.56 

6.59 

6.60 
6.61 


But  we  cannot  determine  all  these  variables  directly  from  hot-wire 
measurements.  A  usual  compromise  is  to  plot  Ae/Aec  against  AeJl/Ae(J  to 
determine  which  variables  may  be  neglected  (ref.  1). 

Since  turbulent  flows  vary  in  both  space  and  time,  useful  insight  into 
turbulent  structure  can  be  gained  from  space-time  correlations  in 
compressible  flows.  Two  approaches  will  be  aescribed.  In  the  first,  let 
us  locate  two  wires  at  points  A  and  B  normal  to  the  flow  direction  such 
that 


eA  =  SAmA  +  S6A0A 

6.62 

*  ®BmB  +  ^6B0B 

663 

Cross-correlating,  we  obtain 
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6.64 

6.65 


This  plots  linearly  in  the  z,  rA  plane  with  intercept 


(  0A0B  “  rB  mB0A^ 


and  slope 

So  that,  if  we  fix  rfi  =  rB1  which  is  a  known  constant,  and  test  with  known 
values  rA  =  rA1  and  rA2  to  obtain  Zj  and  z2,  we  obtain 

^slopeVB1*  const  =  rBimAmB  ‘  mAeB  6  66 

and 

(intercept)  =  con$t  *  0^eB  -  rB1  m^A  6.67 

B  X 

If  we  repeat  procedure  with  a  new  known  constant  rfi2  we  obtain  four 

equations  and  four  unknowns.  To  improve  our  accuracy  we  can  continue 
the  test  sequence  and  determine  the  four  unknowns  graphically  as  shown 
on  the  next  page. 
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wv-v*-. 


intercept 


Then  plot  (intercept,  I)  =  0AOB  rB 


intercept  = 


slope  =  mB6A 


Space-time  correlations  in  compressible  flows  may  be  obtained  with  two 
wires  at  points  A  and  B  in  the  flow  as  follows: 


eA  =  aAmA  +  Pa® a  at  ^me  1 


6.68 


eB  =  aBmB  +  PB0B  at  time  t’ 


6.69 


If  we  adjust  the  anemometer  gains  such  that  rA  =  rB  =  r  and  measure 

(<V  6B>2  ~  <V  gB)2 

<eA*  «8>2  *  <eA-  ®B)2 


we  obtain 

r2  +  r  +  +  ^b  +  c 

Ka  o  =  _ _  - 


'eA®B 


r2  n»i  *  2rA  mAeA  ♦  el 


6.70 


where 


r_  /«A  <*B  r  _  “a 

«/Ba  Ba 


ocA  - 

er  m*m* 


When  the  wires  are  made  of  the  same  material  and  operated  at  the  same 
overheat  ratio  then 


«a  „ 

®A  ~ 


*  c  =  0  &  r\  =  r2 

A 


Furthermore,  as  the  overheat  ->  0,  a/P  -4  0  and  equation  6.70  reduces  to 
the  temperature-temperature  correlation 


0AeB 

eAeB 


At  high  overheat  ratios  a/p  ->1.0,  or  »  1.0,  so  that  if  0  «  m,  or  0  *  m,  the 
equation  reduces  to  the  mass  flux  -  mass  flux  correlation 


mAmB 

mAmB 

To  determined  the  statistical  values  outlined  earlier,  we  have  to  add, 
subtract  and  multiply  the  fluctuating  outputs  and  sometimes  delay  them 
in  time  and  filter  with  frequency.  Addition  and  subtraction  are 
accomplished  by  converting  signal  voltages  to  currents  and  applying 
them  to  resistors  as  shown  in  the  circuit  below. 


In  the  circuit,  A  has  a  high  input  resistance  (»  R.3),  so  that  all  current  goes 
through  R.3,  a  large  voltage  amplification  occurs  at  A  and  the  voltage  at  s 
which  is  equal  to  eQ/A,  tends  to  0.  Now 

e0  =  -  r3  +  ~  )  6  72 

VRl  r2' 

so  that  if  R1  =  R2  then  eQ  =  -(e,  +  e2).  Also,  if  e2  =  0  then  Ri  =  R3  and  eo  = 

-e.  In  these  ways  we  are  able  to  add  and  subtract  hot-wire  voltages.  We 
can  integrate  hot-wire  signals  as  follows 
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c 


s(e=0) 


Again  if  A  ->  es  »  0  and  the  input  impedance  of  A  «>.  For  the 
capacitance  Q  =  Ce0  we  may  write 


dQ  _  ,  _  r  _  ±i 

dt  dt  "  R 


so  that 


and  if  e-  =  a  sin  00  then  eQ  =  -(a/coRC)  cos  cot. 

Differentiation  may  be  accomplished  as  follows 
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i  i 

1 

1  A  _ 

1  ”1 

1 A  ■■ 

s(e=0) 


where  we  may  write 


so  that 


C^-  i  - 
dt  R 


e0-RC^ 


6.73 


6.74 


6.75 


6.76 


In  general,  what  is  termed  the  quarter-square-technique  is  used  for 
multiplication.  In  this  case  two  squaring  circuits  are  used,  one  to  square 
the  quantity  (x  +  y)  and  the  other  to  square  the  quantity  (x  -  y).  The 
outputs  of  the  two  squaring  circuits  are  then  scaled  such  that 


TO 


6.77 


(l/4)[(x  +y)2  -  (x  -  y)2)  =  xy 
The  input  is  usually  scaled  such  that  the  output  is  (l/10)xy.  This  scaling 
is  necessary  to  use  the  full  dynamic  range  of  the  two  input  amplifiers 
without  saturation  of  the  output  amplifier.  Division  can  be  accomplished 
by  placing  a  squaring  network  in  a  feedback  loop  as  shown  below 


The  power  spectrum  of  a  hot-wire  signal  can  be  obtained  by  Fourier 
transformation  of  the  auto-correlation  or  by  spectrum  analysis  using 
inductance-capacitance  band  pass  filters.  However,  filters  with  band 
width  of  1/3  octave,  ie.  21/3  x  center  frequency,  are  not  suitable  for 
measuring  spectra  with  sharp  peaks.  Typical  turbulent  velocity  spectra 
drop  off  rapidly  with  frequency.  But,  in  supersonic  flows,  there  may  still 
be  significant  energy  beyond  50  kHz.  The  sharp  drop  off  can  cause 
changes  in  the  effective  filter  band  width  and,  at  extremely  high 
frequencies,  the  spatial  resolution  of  the  hot-wire  sensing  element 
renders  spectral  measurements  open  to  question.  Although  spectra  of 
higher  order  correlations  have  been  reported  at  low  speeds,  care  must  be 
exercised  to  match  the  phase  shifts  of  the  two  filter  sets,  when  making 
these  measurements  in  high-speed  flows. 


set  ups  are  shown  below 


Intensity  Spectrum 


Cross  Spectrum 


Auto  or  Cross  Correlations 


Details  of  the  properties  of  materials  used  in  the  fabrication  of  hot-wire 
probes  are  given  in  an  attached  list.  Generally  hot-wire  probes  are 
welded  or  soft-soldered  to  the  supports.  Conventional  or  etched 
wollaston  wires  are  used.  Film  gages  have  much  less  directional 
sensitivity  because  of  conduction  to  the  glass  (ie.  their  effective  aspect 
ratio  is  lower).  They  also  have  a  complex  time-constant  so  it  is  difficult  to 
compensate  for  quantitative  measurements.  Some  examples  of  probe 
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bassist 


leads  on  surface 


The  fabrication  sequence  for  double  sensor  with  straight  or  coiled  wires  is 
shown  as  follows: 


Nickel  Support  Pins 


Ceramic  Insert 


1 1 1 


.125  Drill 


U 

i  - 3 

.200  ±  .002 

1  *1 

« — .500 — > 

<—.375—* 

.002 

.004 


.150+  0005 


h -  .675  - H 


Aluminum  or  Stainless  Steel  Body 


4  v 


14.5 


-.10 


Miniature  Co-axial  Cable  Beldon  8700,  28AWG 


(2)  Ni  oxide  coated  pins 

(2)  Ni  pins 

Epoxied  to  insert 
(BB  2101  Epoxy) 


Pin/Ceramic  Assembly 


112 


mwwj 


'  • 


SOLDER 

INSULATE  JOINTS 
WITH  EPOXY 


■e 


14.5  —♦ I 


Cable/Pin/Ceramic  Assembly 


~  .10 

distance  of 
unobstructed 


Preparation  of  Pins  for  Coil  Winding 
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Finally  wind  the  coil,  straighten  bent  pins,  spot  weld  single  wire  to  Ni 
pins,  record  resistances;  check  for  shorts  to  ground,  and  coat  spot-weld 
with  conductive  epoxy  (E3021). 

The  procedure  for  the  construction  of  surface  wire  gages  is  shown  on  the 
following  pages.  These  probes  can  be  used  for  transition  and  separation 
detection  without  disturbing  the  flow.  Surface  flow  direction  can  also  be 
sensed  with  a  pair  of  wires  on  a  single  gage  (ref.  11). 


Stainless  Steel  Tube 
(0.125"  O.D.,  3"  long) 

Shield  of  cable - 


0.125  O.D. 


Electrical  leads 
(0.015"  dia  Ni  Wires) 
Polystyrene  Plug 
(cast  from  resin) 

Core  lead  of  cable 


Epoxy  filling 


Coaxial  Cable 

ending  in  BNC  (female) 

connector. 


Sensor  Wire 

(0.0002"  dia.  Tungsten  Wire) 
spot  welded  to  electrical 
leads  and  solvent  buried 
using  methyl-ethyl-ketone 
as  solvent. 


h-  0.0625 
Surface  Wire  Gage 
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45° 


0.125  O.D. 
-0.116  O.D 


0.08 


Surface  Flow  Direction  Sensor  (ref.  12). 
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Summary  of  Materials 


Section  7  Hot-Wire  and  Laser  Velocimeter  Measurements 


In  the  first  part  of  this  section,  the  results  of  our  hot-wire  freestream 
turbulence  measurements  in  the  M  =  6.0  and  M  =  3.0  high  Reynolds 
number  facilities  will  be  described.  Comparisons  with  previous  data 
obtained  at  much  lower  unit  Reynolds  numbers  in  the  Ames  3.5  foot  HWT 
(fig.  7.1)  are  very  encouraging  as  they  show  superior  flow  quality  in  the 
AFWAL  facility  at  high  unit  Reynolds  numbers.  The  mass  flow 
fluctuations  increase  with  tunnel  total  pressure  and  range  from  0.6  to  1.6 
percent.  Our  previous  Ames  data  were  taken  in  1972  and  after  the 
tunnel  was  converted  to  a  free  jet  test  section.  These  latter  data  were 
taken  in  1974. 

If  we  assume  that  the  disturbances  sensed  by  the  hot-wire  are 
predominantly  sound  waves  radiated  from  the  turbulent  nozzle  wall 
boundary  layers,  then  the  pressure  fluctuation  levels  can  be  estimated 
from  our  hot-wire  data.  Hot-wire  theory  shows  this  assumption  to  be 
consistent  with  linear  mode  diagrams,  (fig.  7.2).  The  results  of  these 
calculations  are  shown  in  fig.  7.3  and  comparisons  made  with  the  Ames 
HWT  and  the  Langley  VDT.  This  figure  clearly  shows  better  flow  quality 
in  the  M  =  6  facility.  But,  sound  is  not  the  only  disturbance  mode,  as 
temperature  spottiness,  probably  due  to  non-uniform  heating  of  the 
supply  gas,  is  not  negligible  (fig.  7.1).  Thus,  the  pressure  level  estimates 
from  the  hot-wire  data  should  be  viewed  as  upper  bounds.  The  actual 
levels  should  be  somewhat  lower.  Direct  pressure  measurements  should 
confirm  this. 


Two  hot-wire  signals  are  shown  in  fig.  7.4;  one  for  a  pressure  of  930  psi, 
the  other  for  a  pressure  of  1860  psi.  These  traces  clearly  show  the 
increased  high  frequency  (smaller  length  scale)  contribution  at  the  higher 
tunnel  total  pressure.  Low-frequency  (large  scale)  contributions  are  also 
apparent  in  both  hot-wire  traces.  In  general,  most  of  the  energy  is 
concentrated  at  low  frequencies.  Auto-correlation  measurements  show 


the  turbulent  integral  length  scales  to  be  of  the  order  of  the  jet  exit 
diameter.  Freestream  turbulence  measurements  have  also  been  made  in 
the  M  =  3  facility.  These  results  are  shown  in  fig.  7.5. 


The  laser  velocimeter  data  were  obtained  across  flat  plate  boundary 
layers  in  the  FDL  M  =  6.0  wind  tunnel.  Zero  pressure  gradient  and 
adverse,  ramp  induced,  pressure  gradient  flows  were  investigated.  In  the 
latter  case,  the  capability  for  both  streamwise  and  vertical  velocity 
measurements  was  demonstrated.  Previously,  problems  with  particle 
arrival  rates,  run  to  run  velocity  variations  and  probe  volume  location 
within  the  boundary  layer  had  been  encountered.  These  problems  had 
led  to  significant  uncertainties  and  scatter  in  the  measured  velocity 
profiles. 


Although  more  costly,  laborious  and  tedious  to  operate,  the  laser 
velocimeter  probably  represents  the  instrument  of  last  resort  for  the 
measurement  of  flow  in  compression  corners,  shock  boundary  layer 
interactions  and  other  large  scale  unsteady  turbulent  flows.  Once  in 
operation,  linear,  non-intrusive  unambiguous  turbulent  velocity  and 
shear  stresses  can  be  obtained  once  seeding  and  sampling  bias  problems 
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have  been  overcome.  Flows  of  most  practical  interest  and  importance 
often  involve  high  turbulence,  flow  separation  and  large  scale 
unsteadiness.  Here  the  inherent  linearity,  non-perturbing,  directionally 
sensitive  properties  of  the  laser  velocimeter  come  to  bear.  However,  it  is 
also  at  these  conditions  that  laser  velocimeter  sampling  bias  comes  to  the 
forefront. 

Most  researchers  agree  that  continuous-wave  mode  signals  are  free  from 
bias.  However,  in  high-speed  wind  tunnel  applications  where  particle 
concentrations  are  low,  individual  realization  processing  is  required.  In 
these  cases  the  potential  errors  attributable  to  sampling  bias  can  become 
significant  at  high  turbulence  levels.  In  a  recent  paper  (ref.  12)  the 
existence  of  "sampling  bias"  in  individual  realization  laser  velocimeter 
measurements  has  been  experimentally  verified  and  shown  to  be 
independent  of  the  sample  rate.  The  results  clearly  demonstrate  that,  for 
the  individual  realization  mode  of  laser  velocimetry,  sampling  bias  exists 
and  that  it  increases  approximately  with  the  square  of  the  turbulence 
intensity.  It  was  also  demonstrated  that  these  bias  effects  are 
independent  of  sampling  rate  provided  the  seeding  concentration  is 
sufficiently  low  to  insure  individual  realization  measurements.  A 
two-dimensional  weighting  of  the  velocity  samples  was  shown  to  be 
effective  in  correcting  the  individual  realization  measurements  for 
sampling  bias. 


With  particle  bias  problems  identified  and  corrected,  the  seeding  problem 
was  the  first  to  be  addressed  during  the  present  work.  A  fluidized  bed  of 


carbon  particles  which  proved  to  be  ineffective  was  replaced  by  an 
atomizer  ahead  of  the  throat.  The  accelerations  and  shear  through  the 
nozzle  were  sufficient  to  break  up  and  provide  an  ample  supply  of 
scattering  centers  in  the  laser  velocimeter  focal  volume.  Signal  to  noise 
was  also  improved  by  the  use  of  pre-amplifiers  on  the  photomultiplier 
tube  outputs.  With  these  improvements,  valid  velocity  data  rates  were 
increased  to  3,000  samples  per  second. 

The  improved  data  acquisition  rate  also  alleviated  the  thermal  growth 
problem  and  its  effect  on  probe  volume  location  within  the  boundary 
layer.  Increased  data  rates  in  the  wall  region  enabled  shorter  run  times, 
reduced  model  heating  and  less  thermal  distortion.  Measurements  can 
now  be  made  through  the  wall  region  without  excessive  model  heating, 
thus  minimizing  measurement  location  uncertainties.  Reliable 
measurements  of  the  axial  velocity  component  were  made  within  0.010 
inches  of  the  wall.  Vertical  velocity  measurements  were  limited  to  0.03 
inches  from  the  wall. 

The  problem  of  run  to  run  edge  velocity  variations  was  first  noticed 
during  preliminary  freestream  shakedown  runs.  These  variations  were 
corrected  by  correlating  each  set  of  shear  layer  measurements  with 
individual  run  total  temperature  measurements.  A  more  elusive  data 
reduction  problem  was  also  identified.  It  concerns  the  relatively  crude 
filter  ranges  on  the  laser  velocimeter  signal  processor.  Since  the  mean 
flow  velocity  and  turbulence  measurements  were  made  with  a  dual  beam 
velocimeter  utilizing  bragg  cell  frequency  shift,  a  stationary  particle 


produced  a  doppler  frequency  of  (f  )  14  mHz  for  example  in  fig.  7.6. 


I 


Thus,  in  the  flow  field,  particles  generated  doppler  frequencies  f  ±  f 

depending  on  their  direction  normal  to  the  moving  fringes.  However, 
even  with  this  offset  errors  can  arise  due  to  incorrect  filter  settings. 
Extra  care  in  their  selection  must  be  exercised  not  only  in  separated  flow 
regions  but  also  in  the  wall  region  of  attached  flows  where  local 
turbulence  levels  are  high.  The  effect  of  16  mHz  and  8  mHz  filter  settings 
are  shown  by  the  solid  and  dashed  lines  on  fig.  7.6.  The  problem  would 
be  more  pronounced  closer  to  the  wall. 
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Fig.  7.7  shows  a  zero  pressure  gradient  flat  plate  velocity  profile.  It  can 
be  seen  that  the  comparison  with  the  Van  Driest  theory  is  excellent.  Fig. 
7.8  shows  a  profile  measured  across  the  ramp  induced  pressure  gradient 

flow  field.  A  direct  comparison  of  ramp  induced  effects  on  the  mean  and 

turbulent  flow  fields  can  be  seen  in  fig.  7.9  where  measurements 
obtained  at  the  same  streamwise  station  are  presented.  These 
measurements  obtained  at  a  station  0.3  6  ahead  of  the  30-deg.  ramp 
clearly  show  retardation  of  the  flow  and  a  significant  increase  in 
turbulence  level  over  a  wide  region.  The  vertical  velocity  profile 
measured  at  the  same  location  ahead  of  the  interaction  is  shown  in  fig. 

7.10.  Local  flow  angularity  profiles  across  the  boundary  layer  have  been 
calculated  from  the  two  component  laser  measurements.  These  results 
(fig.  7.11),  show  that  local  flow  angles  close  to  the  wedge  deflection  angle 
are  present  in  the  shear  layer  just  upstream  of  the  interaction. 

Comparison  of  the  two  turbulence  level  profiles  shows  that  the 


streamwise  turbulent  kinetic  energy  for  the  ramp  flow  is  more  than  three 
times  that  for  the  flat-plate  boundary  layer.  Turbulent  mixing  length 
scales,  calculated  using  local  rms  levels  and  mean  gradients,  are  an  order 
of  magnitude  larger.  Turbulence  levels  based  on  local  mean  flow  values 
exceed  30  percent  in  the  wall  region  so  that  we  can  see  from  section  4 
that  significant  hot-wire  measurement  errors  would  arise.  At  this  high 
intensity,  large-scale  turbulence  results  in  directional  intermittency  of  up 
to  15  percent  ahead  of  the  time-averaged  recirculation  zone.  Hot-wire 
measurement  errors  associated  with  directional  intermittency  are 
discussed  in  detail  in  ref.  7. 

To  conclude,  diagnostic  tools  are  now  available  to  attempt  the 
measurement  of  turbulent  high-speed  flows,  an  area  where 
comprehensive  studies  are  lacking.  However,  measurement  techniques 
must  be  used  with  understanding  and  care  in  appropriate  situations.  All 
too  often  experimental  methods  have  been  stretched  beyond  their 
reliable  limits  and  misleading  results  have  been  reported.  Apparent 
discrepancies  between  measurements  and  calculations  cannot  solely  be 
attributed  to  deficiencies  in  turbulence  modeling  until  reliable 
assessments  have  been  made  of  measurement  errors.  Although  the 
potential  of  laser  velocimetry  for  the  non-intrusive  measurement  of 
mean  velocity,  turbulence  intensity  and  shear  stress  in  complex  wind 
tunnel  flow  fields  has  long  been  recognized,  the  design,  construction  and 
successful  operation  of  systems  in  other  than  small,  closely  controlled 
laboratory  situations  still  presents  a  formidable  engineering  challenge. 
Thus,  it  is  important  that  redundant  hot-wire  and  laser  velocimetry 
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Fig.  7.1  Mass  flow  and  total  temperature  fluctuations 
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Fig.  7.3  Freestream  pressure  fluctuations 
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Fig.  7.7  Flat  plate  axial  velocity  profile 
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List  of  Symbols 


AW 

aW 

C 

Cxy 

C 

E,  e 
e\  Ae 

Aepu  ,AeTt,Ae0 
f 

Gx(t) 

Gxy(t) 

Gr 

I 

i,  Ai 

K 

k2 

k 

L 

M 

m 


Overheat  parameter,  — 

H  2  RwVdl/ 

Rw  '  Rr 

Overheat  parameter,  — — 

Thermal  capacity,  electrical  capacitance 
Modulus  of  cross  spectral  density  function 
Specific  heat 
DC  voltage 
Unsteady  voltages 

Mass  flux,  total  temperature  and  angle  sensitivities 
Frequency 

Power  spectral  density 
Cross  spectral  density 
Grashof  number 
DC  current 
Unsteady  currents 

d  log  Rw 
d  log  Tw 

constant  in  equation  3.13 
Thermal  conductivity 
Wire  length,  turbulent  mixing  length 
Mach  number,  wire  time  constant 
Mass  flow  fluctuation 
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Nu 


P 

P(x),  P(x,y) 
Pr 

q2 

Q(x,y) 

Re 

R 

r 

s 

T,  T 

u,  V,  "w 

U  1,  U2,  u3 
u',  v\  w' 

v,  V 
Zc 

x,  y,  z 


Nusselt  number 
d  log  k 
d  log  Tw 

Pressure 

Probability,  joint  probability  functions 
Prandtl  number 

Defined  inequation  3.3  (turbulent  energy) 
Modulus  of  cross  spectral  density  function 
Reynolds  number 

Wire  resistance,  correlation  function 
Sensitivity  ratio,  fluctuating  wire  resistance 
Sensitivity 
Temperature 

Mean  velocity  components 

Fluctuating  velocity  components  in  section  4 

Fluctuating  velocity  components 

Mean  and  fluctuating  voltage 

Impedence 

Cartesian  Coordinates 


Temperature  coefficients  of  resistance 


l/(,+  V-M2) 

a (Y  -  1)M2 

Ratio  of  specific  heats 

Boundary  layer  thickness 

Displacement  thickness 

d  log  Iw 

Finite  circuit  factor  -  - -  ,  viscous  dissipation 

0  log  Rw 

Recovery  factor 

Temperature  fluctuation,  overheat  parameter  ~ 
Wavelength,  heat  transfer  coefficient 
Viscosity,  wire  time  constant 
Density 

Kinematic  viscosity,  voltage  ratio 

.  Tw"Tr 
Temperature  loading  — — — 

Tr 

Skin  friction 

RMS  value 

Time  delay 

Wire  inclination  angle 

Mean  square  value 

Time  constant 


Subscripts  denote  values  evaluated  at 


c 

f 

g 

m 

n 

o,  t 

r 

w 

oo 

p 

pu 

u 

V 


Compensated 
Film  temperature 
Gas  temperature 
Measured 
Noise 


Tw+Tg 

2 


Total  temperature 
Recovery  temperature 
Wire  temperature 
Freestream 


Due  to  density 
Due  to  mass  flux 


Due  to  axial  velocity 
Due  to  vertical  velocity 


Superscripts  denote 

’  Fluctuating  component,  rms  value 

~  rms  value 
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